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PROJECT  SUMMARY 


The  proposed  research  entitled  'Computational  Modeling  System  for  Deformation  and  Failure  in 
PolycrystaUine  Metals’  will  develop  an  integrated  system  of  experimentally  validated  computational 
models,  for  simulating  the  deformation  and  failure  behavior  of  polyciystalline  metals  and  alloys.  The 
system  will  consist  of  computational  models  and  software  for:  (a)  3D  microstructure  reconstruction  and 
characterization  of  morphological  and  crystallographic  information,  (b)  image-based  microstructural 
computational  models  incorporating  crystal  plasticity  and  fracture  laws,  and  (c)  a  multi-time  scaling  time 
integration  algorithms  for  cyclic  deformation,  to  investigate  the  effect  of  material  microstructure  on 
deformation  behavior  and  fatigue  life.  The  computational  system  is  expected  to  assist  engineers, 
especially  in  the  aircraft  engine  industry,  in  their  product  design  for  enhanced  performance  and  reliability. 

It  is  well  known  that  the  development  of  such  an  integrated  system  is  a  necessary,  yet  challenging 
undertaking.  Conventional  commercial  Finite  Element  Method  (FEM)  codes  have  a  number  of 
shortcomings  when  modeling  real  materials  due  to  lack  of  certain  key  aspects,  such  as  efficiency, 
accuracy  in  the  presence  high  gradients  and  discontinuities,  numerical  stability,  etc.  Important  advances 
are  needed  in  computational  tools  and  methodologies  to  overcome  these  shortcomings.  The  PI  has  gained 
considerable  research  experience  in  the  development  of  advanced  micromechanical  computational  models 
for  material  characterization  and  modeling  that  can  effectively  address  some  of  these  problems  for 
poiycrystalline  and  poly-phase  materials.  An  important  component  of  the  development  of  any 
computational  model  is  its  validation  with  experiments.  Collaborations  of  the  PI  with  Dr.  M.J.  Mills 
(OSU),  Dr.  D.  Dimiduk  (AFRL/MLLLM)  and  Dr.  M.  Uchic  (AFRL/MLLM)  in  his  current  research  [1,2], 
have  provided  this  indispensable  component  to  add  robustness  to  the  models.  The  proposed  research  will 
involve  collaborations  with  Drs.  Mills,  Fraser  of  OSU,  and  Drs.  Dimiduk,  Uchic  and  Parthasarathy  of 
AFRL/MLLM,  for  insight  and  validation  of  the  computational  systems  to  be  developed.  The  research 
partners  have  complimentary  strengths  (hat  will  contribute  towards  the  comprehension  and  overall  goals 
of  this  research.  The  proposed  program  will  complement  ongoing  research  thrusts  in  the  Materials 
Directorate  at  AFRL  and  at  OSU  under  the  Science  and  Technology  Workforce  for  the  21st  Century  (STW~ 
21)  cooperative  agreement  program  between  OSU  and  AFRL. 

Computational  models  will  be  developed  in  the  proposed  program  to  incorporate  accurate 
microstructural  morphology  and  crystallographic  information  for  analyzing  structure-sensitive 
deformation  and  fatigue  behavior.  Primary  focus  will  be  on  aerospace  metals  and  alloys,  e.g.  IN-1  DO,  Ti- 
6242  and  Ti-6Al-4V,  with  a  motivation  to  improve  performance  and  reliability  of  aerospace  engine 
materials.  A  comprehensive  approach  will  be  pursued,  coupling  the  following  elements: 

(i)  Reconstruction  of  2D  and  3D  microstructural  models  with  real  morphological  features  and 
crystallographic  orientations  from  quantitative  metallography  combining  focused  ion  beam  (FIB) 
and  electron  back-scatter  detector  (EBSD)  in  orientation  imaging  microscopy,  microstructural 
characterization,  and  statistical  methods. 

(ii)  Development  of  experimentally  validated  finite  deformation  rate-dependent  crystal  plasticity 
models  with  size  effects,  for  analyzing  deformation  and  creep  in  polycrystalline  aggregates. 

(Hi)  Development  of  an  experimentally  validated  cohesive  zone  model  for  inter  and  intra-granular 
fracture  initiation  and  growth  for  fatigue  crack  modeling. 

(iv)  Development  of  a  stable,  accurate  and  efficient  image-based  micromechanical  finite  element 
model,  for  crystal  plasticity  and  damage,  incorporating  real  morphological  and  crystallographic 
features.  A  novelty  of  this  model  would  be  in  the  powerful  adaptive  techniques  for  automatically 
changing  the  element  resolution  and  topology  with  evolving  strain  localization  and  damage. 

(v)  Development  of  multi-scaling  algorithms  in  the  time  domain  for  compression  and  localization 
in  cyclic  deformation  analysis. 

In  summary,  the  proposed  research  will  develop  experimentally  validated  innovative  high  fidelity 
computational  tools  to  provide  improved  predictive  capabilities  for  fatigue  failure  analysis  in  materials 
with  aerospace  applications.  The  program  is  expected  to  advance  the  state  of  the  art  in  analysis 
methodology  for  high  confidence  in  component  performance  and  life  calculations. 
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1. 


STATEMENT  OF  OBJECTIVES 


Following  the  Air  Force’s  thrusts  in  advanced  military,  aerospace  and  space  exploration  systems,  the 
proposed  research  is  aimed  at  developing  an  integrated  system  of  computational  models  and  predictive 
tools  that  can  be  utilized  to  enhancing  performance  and  reliability  of  polycrystalline  metals  and  alloys 
used  in  aerospace  applications.  These  materials,  e.g.  IN- 100  or  Ti-6242  and  Ti-6A1-4V  are  used  in  aircraft 
engine  components.  Understanding  the  role  of  material  microstructure,  at  the  length  scale  of  grains  and 
polycrystalline  aggregates,  on  the  deformation  and  failure  characteristics  of  the  material  is  critical  to  the 
reliable  design  of  components  for  high  performance  applications.  Such  an  understanding  requires  an 
analysis  framework  that  can  predict  inhomogeneities  in  time-dependent  plastic  flow  under  fatigue  and 
creep  conditions.  Naturally,  that  sets  a  requirement  for  representing  the  real  microstructure  and  defects, 
within  the  plasticity  analysis  tools.  A  robust  design  methodology  must  also  link  variabilities  involved  at 
all  length  scales  that  can  affect  the  components  in  service  performance. 

Developing  advanced  computational  tools  to  unravel  these  complex  relations  will  be  the  objective  of  the 
proposed  research.  Most  of  the  current  efforts  in  this  area  employ  commercial  FEM  codes  like  ABAQUS 
[3],  Many  of  these  implementations  suffer  from  severe  limitations  in  efficiency,  accuracy  and  resolution 
in  regions  of  high  gradients,  as  well  as  numerical  stability,  for  modeling  real  materials.  The  PI  has  been 
involved  in  the  development  of  advanced  micromechanical  computational  models  for  material 
characterization  and  analysis  of  deformation  and  failure  in  poly-phase  and  polycrystalline  materials.  This 
experience  will  be  utilized  in  the  proposed  research  to  develop  numerically  stable  and  robust,  image- 
based  micromechanica!  computational  models  of  plasticity  and  damage,  incorporating  real  morphological 
and  crystallographic  features.  The  goal  will  be  to  investigate/predict  the  dependence  of  cyclic  deformation 
behavior  and  fatigue  crack  evolution,  on  the  material  microstructure.  The  model  development  will  be 
intricately  coupled  with  experimental  validation,  for  which  the  material  system  will  concentrate  on 
Titanium  and  Nickel  alloys  with  hep  and  fee  crystalline  structure  respectively.  The  program  will 
complement  ongoing  research  thrusts  on  materials  development  and  structural  design  at  the  Materials 
Directorate  of  AFRL/MLLM  under  the  leadership  of  Dr.  D.  Dimiduk.  Ongoing  collaborations  with 
Professors  M.J.  Mills  and  H.  L.  Fraser  atOSU  and  Drs.  D.  Dimiduk,  T.A.  Parthasarathy  and  M.  Uchic  at 
AFRL/MLLM  will  be  augmented  under  the  auspices  of  this  research  program.  A  comprehensive 
approach  will  entail  the  following  activities. 

i.  Construction  of  2D  and  3D  Microstructural  Models:  This  task  will  create  computer  models  of  the 
material  microstruelure,  incorporating  morphological  features  like  grain  geometry  and  neighbor 
description,  and  crystallographic  orientations.  Quantitative  metallography,  combining  focused  ion 
beam  (FIB)  techniques  for  material  sectioning  and  electron  back-scatter  detector  (EBSD)  in 
orientation  imaging  microscopy  (OIM),  will  provide  the  necessary  microstructural  information  and 
data.  Such  experimental  activity  is  already  in  progress  at  OSU  and  AFRL  facilities  under  the 
existing  Science  and  Technology  Workforce  for  the  2]*'  Century  (STW-2J)  cooperative  agreement 
between  AFRL  and  OSU.  Subsequently,  microstructural  characterization,  statistical  and  stereological 
methods  and  computational  algorithms  will  be  explored  for  effective  3D  construction. 

ii.  Development  of  Experimentally  Validated  Constitutive  and  Damage  Models:  In  this  task,  finite 
deformation  rate-dependent  crystal  plasticity  models  with  size  effects  will  be  developed,  for 
analyzing  cyclic  deformation  and  creep  in  polycrystalline  aggregates.  In  addition,  special  cohesive 
zone  models  with  traction-displacement  laws  will  be  developed  for  inter  and  intra-granular  fracture 
initiation  and  growth  during  cyclic  deformation  and  fatigue  cracking.  The  models  will  be  validated 
with  experiments  for  plastic  deformation  and  creep  in  polycrystalline  aggregates  and  evolution  of 
fatigue  fracture.  This  research  will  be  integrated  with  ongoing  experimental  research  programs  at 
OSU  (Mills,  Fraser)  and  AFRL/MLLM  (Dimiduk,  Parthasarathy,  Uchic)  on  state  of  the  art 
mechanical  testing  and  material  characterization. 


4 


iti*  Development  of  Adaptive  Image  based  Micromechanical  Computational  Models;  This  task  will  build 
novel  finite  elements  models  that  evolve  from  the  material  microstructure  and  have  special 
characteristics  of  numerical  stability,  accuracy  and  resolution  and  efficiency,  for  highly  non- 
homogcneoiis  fields*  In  particular,  a  grain-based  finite  element  model  (GBFEM)  manifesting  the 
material  crystalline  structure  and  morphology  will  be  developed  using  computational 
micro  mechanics  for  finite  deformation  rate-dependent  crystal  plasticity  and  micro  structural  damage 
evolution.  Adaptivity  of  the  computational  domain  is  a  necessary  feature  to  enhance  efficiency, 
especially  in  the  presence  of  evolving  localized  deformation  and  cracking.  GBFEM  will  be  endowed 
with  adaptive  capabilities  to  identify  and  impart  preferential  resolutions  to  these  regions  by  local 
enrichment*  It  will  also  admit  changing  topology  due  to  the  onset  and  growth  of  microcracks  through 
adaptive  introduction  of  cohesive  zone  models  of  fracture.  Time  stepping  algorithms,  currently  used 
in  commercial  codes  for  the  simulation  of  cyclic  deformation  during  fatigue,  have  serious  limitations 
with  respect  to  the  number  of  cycles  that  can  be  simulated.  An  Eidaptive  multi- time  scaling 
methodology  is  currently  being  developed  by  the  PI  to  overcome  this  limitation  through 
homogenization  and  localization  of  computations  in  the  time  domain  [4].  This  multi- time  scaling 
methodology  will  be  incorporated  in  the  computational  modeling  framework.  Finally,  to  facilitate 
large  scale  computing  of  real  materials,  the  models  will  be  coupled  with  high  performance 
computing  tools,  incorporating  parallel  structure  with  dynamic  load  balancing  strategies* 

In  summary,  this  research  will  introduce  new  directions  and  advance  the  state  of  the  art  in  computational 
modeling  methodologies  that  can  provide  an  improved  level  of  understanding  of  the  factors  that  cause 
fatigue  failure  in  metallic  materials* 


2*  MOTIVATION  AND  SIGNIFICANCE 

Advanced  military  aerospace  systems  like  the  Joint  Shake  Fighter,  unmanned  missile  technologies,  future 
military  aircrafts  as  well  as  improved  space  access  systems,  are  thrust  areas  of  priority  at  the  Air  Force 
[4].  Major  goals  of  these  programs  are  to  drastically  reduce  hardware  and  support  costs,  improve 
reliability,  expand  reusability  and  to  enhance  thrust- to- weight  performance.  Improved  materials  having 
high  specific  properties,  e.g.  low  weight,  high  stiffness,  fatigue  strength,  and  fracture  toughness  at  a  range 
of  service  temperatures  are  essential  for  achieving  these  lofty  structural-efficiency  goals*  The  need  for  the 
proposed  research  is  derived  from  these  thrusts  to  improve  space  and  aero-propulsion  systems* 

Progressive  developments  in  advanced  materials  have  contributed  tremendously  to  the  enhancements  in 
aerospace  applications,  especially  in  the  thrust- to- weight  ratio  of  aero-gas  turbine  engines*  This  has 
mainly  been  possible  through  the  introduction  of  titanium  and  nickel  based  alloys  like  super-alloy  IN- 100 
and  Ti-6242  and  Ti-6A1-4V,  which  have  provided  enabling  technology  for  enhanced  performance  and 
reliability.  The  design  and  development  of  lightweight,  high  temperature-creep  resistant  aerospace  engine 
components  like  fan  blades,  high-pressure  compressors  and  high-pressure  turbines  (figure  la)  have 
benefited  significantly  from  these  materials.  Their  contribution  is  vital,  not  only  in  terms  of  desirable 
properties,  but  also  improved  reliability*  Continuing  growth  of  the  aerospace  industry  demands  further 
research  and  development  of  advanced  materials  having  heterogeneous  or  graded  microstructures,  such  as 
the  dual -microstructure  turbine  disks.  The  proposed  research  is  motivated  by  the  need  for  rational  design 
procedures  to  select  ' optimal  material  microstructures  ’  that  wifi  enhance  the  material’s  mechanical 
behavior  and  failure  characteristics* 

Aerospace  components  tire  exposed  to  cyclic  loading  conditions  due  to  start  up  and  shut  down  processes 
or  load  reversals  during  service.  In  many  cases,  this  results  in  their  fatigue  or  time  dependent  fracture* 
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Fatigue  failure  in  the  microstructure  primarily  evolves  in  three  stages,  viz.  (i)  crack  nucleation  due  to 
inhomogeneous  plastic  flow  or  grain  boundary  failure,  (ii)  subsequent  crack  growth  of  difficult  to  predict 
microstructurally  small  cracks,  by  cyclic  stresses  at  a  rate  dictated  by  plasticity,  stress  range  and 
environment,  and  (iii)  coalescence  of  cracks  into  the  faster  crack  propagation  regime.  The  phenomena  of 
high-cycle  and  low-cycle  fatigue  have  been  traditionally  characterized  using  macroscopic  parameters  like 
applied  stresses,  cyclic  frequency,  loading  waveform,  hold  time  etc.,  as  well  as  statistical  distributions  of 
fatigue  life  and  fatigue  strength  [6-1 1].  Fatigue  design  by  total  life  approaches  includes  the  stress-life  or 
S-N  approach  where  the  stress  amplitude  versus  life  is  determined,  and  the  strain-life  approach  e.g.  the 
Coffin-Manson  rule,  where  the  number  of  cycles  to  failure  is  determined  as  a  function  of  plastic  strain. 
Predicting  the  time  to  initiate  a  fatigue  crack  and  propagate  it  to  failure  has  been  the  prime  objective  of 
these  models.  In  the  stress-life  approach,  the  applied  stresses  are  nominally  elastic  and  the  number  of 
cycles  to  failure  is  large  as  in  high  cycle  fatigue,  while  in  the  strain-life  approach  components  undergo 
significant  plastic  straining  and  crack  propagation.  The  total  life  approaches  have  been  adjusted  for  notch 
effects  using  fatigue  strength  reduction  and  for  variable  amplitude  fatigue  e.g.  in  the  Pal mgren- Miner  rule 
of  cumulative  damage.  In  all  cases,  microstructural  effects  are  represented  only  by  shifts  in  such  data 
curves  after  extensive  testing. 

Alternatively,  the  defect  or  damage  tolerant  approaches  determine  fatigue  life  as  the  number  of  cycles  to 
propagate  a  crack  from  a  certain  initial  size  to  a  critical  size.  These  are  determined  from  threshold  stress 
intensity,  fracture  toughness,  limit  load,  allowable  strain  or  allowable  compliance.  The  models  assume  the 
presence  of  a  flaw  in  the  structure  and  predict  life,  using  laws  like  the  Paris  law  [12].  Fatigue  crack 
advance  has  been  conventionally  based  on  linear  elastic  fracture  mechanics  analysis  related  to  the 
concepts  of  similitude,  in  which,  stress  intensity  factors  uniquely  characterize  fatigue  crack  growth.  A 
number  of  studies  on  fatigue  behavior  of  commercial  alloys  [13,14]  have  demonstrated  the  influence  of 
deformation,  damage  mechanisms,  creep,  oxidation  and  microstructural  instabilities  on  cyclic  life.  Creep- 
fatigue  life  predictions  models  have  been  developed  where  the  interaction  between  creep  and  fatigue  has 
been  accounted  for  by  a  frequency  factor  [15],  strain  range  partitioning  [16],  and  damage  functions  based 
on  hysteresis  energy  [17].  Although  these  models  have  worked  well  for  specific  alloys  under  specific  test 
conditions,  a  lack  of  underlying  physics  and  microstructure  based  considerations  plague  their  portability 
to  other  materials  and  load  conditions. 

The  advanced  materials  used  in  many  aerospace  components  have  complex  multi-colony,  multi-phase 
polycrystalline  aggregates  (figure  lb).  The  mechanical  behavior  and  fatigue  failure  response  are 
intricately  governed  by  microstructural  features  that  include  morphological  and  crystallographic 
characteristics,  e.g.  shape,  size  and  location  of  phases  in  the  colony  structure,  relative  colony  size  and 
locations,  crystal  orientations  and  misorientations,  grain  boundary  geometry  etc  (figure  lc).  These 
intrinsic  attributes  alone  would  be  complex  enough,  however  even  further  complexity  derives  from  these 
features  interacting  with  extrinsic  defects  and  design  features.  In  the  recent  years,  detailed 
micro  mechanical  models  have  been  used  to  understand  damage  mechanisms  leading  to  fatigue  failure  and 
to  throw  light  on  the  stochastic  nature  of  fatigue  phenomena.  However,  morphological  and 
crystallographic  heterogeneities  yielding  anisotropic  properties  and  localized  non-homogeneous 
deformation  impose  severe  challenges  to  these  computational  models.  While  some  computational  models 
of  heterogeneous  materials  have  predicted  the  stress-strain  behavior  with  reasonable  accuracy  [13,14] 
their  competence  in  predicting  microstructural  variability  and  effects  on  ductility  and  fatigue  failure  is  far 
from  mature.  Crystal  plasticity  theories  with  explicit  grain  structures  are  effective  in  predicting  localized 
cyclic  plastic  strains  [15].  But  their  use  in  accurate  damage  prediction  by  inter-  and  trans-granular  and 
grain  boundary  decohesion  is  a  difficult  task,  due  to  interactions  of  grain  orientations,  micro-damage, 
strain  localization  and  other  factors.  Studies  on  deformation  modeling  [16]  have  established  that  preferred 
grain  orientations  due  to  crystal  rotation  facilitate  slip  transmittal  across  grain  boundaries  and  are  a 
dominant  cause  for  material  flow  anisotropy  leading  to  failure.  Metallographic  analyses  indicate  that 
aside  from  texture,  highly  non-homogeneous  deformation  at  inter-  and  intra-granular  levels  can 
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significantly  affect  the  macroscopic  softening  and  failure  process  [17],  Experimental  studies  [18]  suggest 
that  the  growth  of  crystallographic  micro-slip  bands  along  active  slip  systems  of  plastic  flow  is  a  likely 
mechanism  for  localized  instability.  These  micro-bands  originate  from  plastic  How  concentration  along 
preferred  slip  directions  and  from  compatibility  requirements  between  interacting  grains.  They  continue 
to  grow  across  grain  boundaries  due  to  grain  structure  instability  and  eventually  manifest  as  macroscopic 
shear  bands.  The  interaction  of  microscopic  shear  bands  with  transverse  grain  boundaries  also  leads  to 
grain  boundary  micro-cracking,  which  grow  in  size  and  merge  to  cause  fracture. 


Turbines 


Figure  I:  (a)  A  gas  turbine  engine t  (h)  optical  micrograph  oc/0  forged  Ti-6242  alloy  used  in  the 
compressor  rotor \  (c)  fatigue  crack  growth  in  01 M  image. 


The  difficulties  are  further  compounded  by  the  diversity  of  mechanical  testing  conditions.  While  simple 
testing  conditions,  e.g.  uniaxial  tension-compression  tests  at  constant  load  amplitudes  [19,20,21]  have 
helped  in  identifying  mechanisms  of  crack  evolution  processes,  their  unambiguous  extrapolation  to 
various  materials  is  questionable* 


In  facing  up  to  the  challenges  faced  in  the  accurate  and  efficient  analysis  and  design  of  advanced 
materials  in  aerospace  components,  the  proposed  research  will  launch  a  research  program  to  develop  an 
experimentally- validated  computational  simulator  for  predicting  material  behavior  and  failure 
characteristics.  Micros tructural  image -based  adaptive  computational  tools  will  be  developed  for  a  better 
understanding  of  the  role  of  material  microstruoture  on  mechanical  behavior.  The  models  will  involve 
poly  crystalline  plasticity  for  multiple  phases  and  crack  growth  laws.  Multi -time  scaling  facilities  will  be 
developed  to  enable  simulation  of  large  number  of  cycles  to  failure.  This  research  is  expected  to  provide 
novel  tools  and  superior  technology  for  better  utilization  of  advanced  materials  in  Air  Force  and  space 
applications. 
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3* 


BACKGROUND  &  PRIOR  RESEARCH  OF  THE  PI  IN  RELATED  AREAS 


In  his  past  and  current  research,  the  PI  has  been  involved  in  development  of  image-based  multiple  scale 
computational  models  for  mechanical  behavior  and  failure  analysis  of  heterogeneous  poty-phase  and 
polycrystalline  materials,  and  for  their  micros  true  tural  characterization.  Some  of  this  research  in  related 
fields  is  summarized  below. 


3.  /  Quantitative  Metallography  and  Material  Characterization  of  Polyphase  Microstructure 
3D  microstructures  were  constructed  for  characterization  and  modeling  by  serial  sectioning  of  particle- 
reinforced  composite  (AfSiC  and  Al-Si  systems)  microstructures  in  [22,23,24].  The  3-D  microstructure 
models  were  computationally  constructed  by  sequentially  stacking  a  series  of  scanning  electron/optical 
micrographs  (figure  2). 


Serially  Stacked 
Micrographs 


Computer 
Visualization  amt 

Integration 


Computer  Simulation 
of  Microstructure  & 
Tessellation 


Figure  2:  Microstructure  for  SiC  particle  reinforced  aluminum  alloy :  (a)  a  Parallel  section  micrographs 
obtained  by  serial  sectioning;  (b)  computer  image  created  by  serially  stacking  micrographs;  (c)  simulated 
equivalent  microstructure  tessellated  into  a  Voronoi  mesh , 

The  microstructures  were  tessellated  into  a  network  of  2-D  and  3-D  Voronoi  cells  (figure  2c).  Methods  of 
image  analysis  have  been  developed  using  spatial  analysis  techniques  such  as  tessellations,  and  statistical 
functions  like  spatial  correlation,  intensity,  and  distribution  functions.  Microstructures  were  characterized 
to  investigate  the  role  of  morphology  on  the  initiation  and  propagation  of  damage.  Statistical  analyses 
were  conducted  with  marked  correlation  functions  to  establish  optimal  sizes  of  the  microstructural 
representative  volume  elements  or  RVE’s. 


3.2  Microstructu  ral  Characterization  of  Polycrystalline  Microstructures  using  FIB/E  BSD  Systems 
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In  this  recent  work  on  the  STW21  project  [2],  a  new  technique  utilizing  an  automated 
approach  of  combining  a  Focused  Ion  Beam  (FIB)  and  Electron  Back -Scatter  Detector  (EBSD) 
has  been  implemented  for  the  collection  of  3-D  orientation  data  that  can  be  used  to  accurately  model  grain 
and  sub-grain  particle  structures.  The  process  is  capable  of  sectioning  a  sample  with  inter-slice 
thicknesses  as  low  as  50  nanometers  and  capturing  an  orientation  map  of  each  slice  (see  figure  3a), 
Reconstruction  of  the  3-D  orientation  data  involves  combining  or  “stacking”  the  2-D  slices  obtained  by 
the  FIB -EBSD  process  (figures  3b  and  c).  A  special  interactive  program  called  Micro- Mesher  is 
developed  in  this  work,  which  uses  the  orientation  data  and  information  to  define  microstructural  features 
such  as  grains  and  second  phase  particles-  Grain  boundaries  are  typically  tortuous  and  have  local  position 
and  orientation  fluctuations.  Micro- Mesher  approximates  these  spatially  complex  boundaries  with  line 
segments  by  using  an  error-per-unit- length  algorithm  as  set  by  the  user  (figure  3b),  Important 
microstructural  statistics  like  size,  number  of  neighboring  grains,  orientations  and  misorientations,  second 
phase  particles  size  and  inter  particle  spacing,  that  are  used  to  define  and  characterize  the  3-D 
microstructure  are  also  calculated  by  Micro- Mesher.  The  3-D  statistical  information  gained  from  this 
process  improves  the  ability  to  accurately  characterize  the  microstructure  [2]. 


Figure  3:  (a)  A  focused  ion  beam  grain  image  of  a  Ti-6242  microstructure  from  an  OIM  scan  (b)  Grain 
map  after  approximation  by  Micro- Mesher  and  sequential  stacking,  and  (c)  a  3D  computer  image  created 
by  serially  stacking  micrographs > 


3.3  The  Voronoi  Cell  Finite  Element  Model  (VCFEM)  for  Micromechanical  Modeling 
Accurate  micro-mechanical  modeling  of  failure  in  complex  heterogeneous  microstructures  requires  very 
high-resolution  models.  A  powerful  Voronoi  Cell  Finite  Element  Method  (VCFEM)  has  been  developed 
in  [25-30]  for  accurate  image-based  modeling  of  non-uniform  heterogeneous  microstructures. 
Morphological  arbitrariness  in  dispersion,  shape  and  size  of  heterogeneities,  as  acquired  from  actual 
micrographs  (figure  2c  and  3a),  are  readily  modeled  by  this  method.  The  VCFE  model  naturally  evolves 
by  tessellation  of  the  microstructure  into  a  network  of  multi-sided  Voronoi  polygons  (figure  4a).  Each 
Voronoi  cell  with  embedded  heterogeneities  (particle,  fiber,  void,  crack  etc.)  represents  the  region  of 
contiguity  for  the  heterogeneity,  and  is  treated  as  an  element  in  VCFEM.  VCFEM  elements  can  be 
considerably  larger  than  conventional  FEM  elements  and  incorporate  a  special  hybrid  FEM  formulation. 
Incorporation  of  known  functional  forms  from  analytical  micro  mechanics  substantially  enhances  its 
convergence.  A  high  level  of  accuracy  and  efficiency  has  been  achieved  with  VCFEM.  It  can  also  provide 
an  excellent  link  between  quantitative  metallography  and  mechanical  response  [28,29].  Successful 
applications  of  2D-small  deformation  VCFEM  have  been  made  in  thermo-elastic  problems  [27],  elastic- 
plastic  problems  of  composite  and  porous  materials  [23,25]  and  damage  initiation  by  particle  cracking 
and  interfacial  debonding  [28,29,30].  An  adaptive  VCFEM  has  been  developed  in  [25]  where  optimal 
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improvement  is  achieved  by  adaptive  enrichment  of  the  displacement  and  stress  fields.  In  figure  3c,  the 
stress  in  a  fiber- reinforced  composite  microstructure  by  the  VCFEM  model  (figure  4a)  is  compared  with 
those  obtained  by  a  high  resolution  FEM  (figure  4b).  Grain  level  anisotropy  and  texture  and  accurate 
representations  of  the  grain  structure  are  still  unresolved  issues  with  this  method  that  will  be  addressed  in 
the  proposed  effort. 


(a) 


Figure  4:  (a)  VCFEM  mesh;  (h)  a  conventional  FEM  (AN SYS)  mesh  for  a  composite  microstructure  with 
dispersed  inclusions;  (c)  Microscopic  stress  variation  by  VCFEM  and  a  high  resolution  ANSYS  mesh. 


3*4  The  Voronin  Celt  Finite  Element  Model  for  Microstructural  Damage  Modeling  by  Interfacial 
Debonding,  Particle  Fragmentation  and  Matrix  Cracking  in  Composite  Microstructures 

The  initiation  and  growth  of  damage  by  interfacial  decohesion  in  fiber-polymer  matrix  composites  is 
analyzed  in  [30]*  The  VCFEM  model  is  augmented  with  cohesive  zone  constitutive  relations  to  model  the 
progress  of  tnlerfactal  debonding.  Inter  facial  traction  increases  with  separation  to  a  maximum  and 
subsequently  subsides  to  zero  traction,  signaling  debonding  in  cohesive  zone  models.  A  combined 
experimental-computational  study  is  conducted  with  polymer  matrix  composite  specimens  with  excellent 
agreement  is  obtained  between  experiments  and  the  simulations.  Numerical  simulations  of  various 
multiple- fibers  microstructures  with  different  morphologies  (dispersions,  volume  fractions  and  size 
distributions)  are  conducted  to  understand  their  influence  on  the  decohesion  process.  Figure  5a  shows  a 
typical  nonuniform  microslructure  of  a  fiber  reinforced  composite  that  is  tessellated  into  a  Voronoi  ceil 
mesh.  Figure  5b  shows  the  corresponding  debonding  pattern  at  the  fiber- matrix  interface.  The  debond- 
induced  damage  is  found  to  be  very  sensitive  to  the  local  morphology  e.g.  inter-fiber  spacing.  It  is  shown 
that  the  initiation  and  especially  propagation  of  debonding  depends  not  only  on  the  total  cohesive  energy, 
but  also  on  the  shape  of  the  traction-displacement  curve.  The  model  is  used  to  study  the  influence  of 
various  local  morphological  parameters  on  damage  evolution  by  interfacial  debonding. 

The  Voronoi  cell  FEM  has  also  been  developed  to  adaptively  simulate  microstructural  damage 
induced  by  particle  fragmentation  and  subsequent  ductile  fracture  of  the  matrix.  An  experimental- 
computational  study  was  conducted  to  understand  damage  phenomenon  that  causes  failure  in  SiC- 
reinforced  aluminum  alloys  [28,29,31].  Through  a  blend  of  mechanical  testing,  microstructural 
characterization  and  micro-mechanical  analysis,  quantitative  aspects  of  microstructural  morphology  that 
are  critical  to  damage  nucleation  was  investigated.  Interrupted  testing  prior  to  complete  failure  of  a  tensile 
specimen  was  conducted  to  gage  the  influence  of  microstructure  on  the  dominant  damage  path.  Various 
characterization  functions  were  generated  and  sensitivity  analysis  was  conducted  to  explore  the  influence 
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of  morphology  on  damage.  VCFEM  was  used  to  model  damage  initiation  by  particle  cracking  and 
propagation  of  damage  through  multiple  particle  fragmentation  into  a  dominant  crack.  VCFEM 
simulations  compared  excellently  with  experimental  observations  as  shown  in  figure  5c  and  5d. 

The  recent  extension  of  VCFEM  is  in  the  coupling  of  damage  modes  by  particle  cracking  and  ductile 
rupture  of  the  matrix.  A  void  initiation  growth-based  continuum  plasticity  model  (Gurson-Tvergaard 
model  for  porous  plasticity)  is  used  to  manifest  ductile  fracture  in  the  matrix.  The  VCFEM  model 
adaptively  evolves  in  (his  augmentation  to  selectively  provide  high  resolution  in  regions  of  plastic  strain 
localization  and  high  void  concentration.  Figure  5e  shows  the  contour  plots  showing  localized  plastic 
strain  and  ductile  fracture  in  a  simulation  with  particle  and  matrix  cracking. 


(a)  (*>) 


(c) 


(d) 


(e) 


Figure  5:  (a)  A  tessellated  mesh  of  Voronoi  cell  elements  for  a  nonuniform  fiber  reinforced  polymer 
matrix  composite  microslructure,  (h)  A  localized  damage  path  induced  by  interfacial  debonding,  (c) 
Experimental  micrograph  of  an  Al-SiC  showing  cracked  particles,  (d)  VCFEM  simulation  of  the 
equivalent  micrograph  showing  damage  and  equivalent  plastic  strain,  (e)  Contour  plot  of  plastic  strain  in 
a  VCFEM  microslructural  simulation  with  particle  and  matrix  cracking. 


3.5  Multi-Scale  Simulations  with  the  Micromechanical  Models 

Composite  structures  inherently  consist  of  multiple-scales,  e.g.  macro,  meso  and  micro-scales.  In  this 
effort,  an  adaptive  multi-level  computational  system  has  developed  for  adaptively  creating  a  hierarchy  of 
computational  sub-doniains  that  provides  necessary  resolution  at  different  scales  [32-37].  Macroscopic 
analysis  is  done  by  conventional  finite  element  method,  while  the  Voronoi  cell  finite  element  model  is 
used  for  meso-  and  microscopic  analyses.  Error  indicators  are  formulated  to  regulate  element  refinement 
and  inter-scale  transition  from  macroscopic  regions  to  microscopic  "hotspots"  of  stress  or  strain 
localization.  The  ultimate  goal  is  to  have  an  automated  computational  model  that  will  model  failure  as  a 
multi-scale  phenomenon  of  damage  initiation  and  propagation.  The  adaptive  multi-level  model  has  been 
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conducted  for  multi-scale  problems  in  elasticity  problems  without  damage  in  [32,33]  and  in  structures 
with  evolving  damage  by  interfacial  debonding  in  the  microstructure  [38,39]. 


Homogenize 


Localize 
Across  Scales 


Interfacial  debonding  in  the 
Microstructure 


Microscale  analysis  with 
VCFEM 


Composite  Structure 


Microscale  RVE 


Figure  6:  Multi-scale  modeling  of  debonding  based  failure  in  a  fiber-reinforced  bonded  joint  with  a  two- 
way  coupled  multi-level  computational  model 


3.6  Modeling  Creep  and  Cyclic  Deformation  in  Polycrystalline  Ti  Alloys  and  US  LA  Steels 

The  PI  has  been  recently  involved  in  two  parallel  experimental -computational  studies  on  analysis  of  creep 
deformation  in  polycryslaliine  Ti-6242  alloys  [1]  and  cyclic  deformation  in  HSLA  steels  [40,41],  in  these 
studies j  a  time-dependent  crystal  plasticity  model  is  incorporated  in  ABAQUS  UMAT  [3],  to  account  for 
large  strains,  material  anisotropy  and  time-dependent  nature  of  plasticity,  integral  to  creep.  The  finite 
element  models  incorporate  crystallographic  orientation  distributions  directly  from  OIM  data.  Statistically 
equivalent  orientation  distribution  is  assigned  to  the  finite  element  model  through  a  special  orientation 
assignment  method  from  OIM  images.  In  this  method,  sets  of  Euler  angles  are  generated  according  to  area 
fraction  of  orientation  distribution  function  (ODF)  values  in  the  OIM,  and  used  to  generate  the  orientation 
distributions  in  the  computational  model  as  shown  in  figure  7.  A  genetic  algorithm  (GA)  based  multi¬ 
variable  optimization  method  is  performed  to  calibrate  crystal  plasticity  parameters  for  various  slip 
systems,  from  a  limited  set  of  experimental  results. 

The  study  on  polycrystals  of  a-Ti-6Al  and  Ti-6242  alloys  [1]  is  aimed  at  gaining  an  understanding  of 

the  role  of  material  microstructure  on  dwell  behavior,  through  a  combination  of  modeling  and 
experiments.  The  inelastic  deformation  by  slip  in  hexagonal  closed  packed  (hep)  Ti-alloys  is  highly 
anisotropic  due  to  difference  in  deformation  resistances  in  different  slip  systems.  An  important  outcome 
of  the  study  is  the  understanding  of  load  shedding  and  local  stress  rising  phenomenon  in  polycrystalline 
aggregate  of  Ti  alloys.  Creep  and  dwell  cyclic  loading  can  lead  to  local  crack  initiation  at  critical 
locations  within  the  microstructure  with  large  orientation  mismatch.  It  is  observed  that  the  stress 
concentration  is  significantly  affected  by  the  material  rate  sensitivity  and  the  peaks  increase  considerably 
with  time.  Considerable  local  yielding  leading  to  plastic  strain  is  realized  in  the  microstructurc,  even  at 
macroscopic  stress  levels  lower  than  the  yield  stress.  Significant  local  stress  concentrations  evolve  with 
time  in  the  material  microstructure  due  to  compatibility  constraint  at  the  interface. 
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Figure  7:  (a)  Micrograph  of  a  forged  a/p  Ti-6242  alloy ,  (h)  contour  plot  of  the  experimental  (0001)  pole 
figure ;  (c)  transformed  experimental  pole  figure  from  contour  plot  to  points  in  the  discretized  regions  of 
the  projected  plane ;  (d)  statistically  equivalent  orientations  obtained  from  the  orientation  assignment 
method  (e)  orientation  distribution  for  the  FEM  simulations  of  polycrystalline  TF6AL 


For  the  understanding  the  behavior  of  HSLA  stee!sr  a  study  combining  experiments  with  finite  dement 
based  simulations  is  conducted  to  develop  a  crystal  plasticity  model  for  prediction  of  the  cyclic 
deformation  behavior  [40,41].  The  experiments  involve  orientation  imaging  microscopy  (OIM)  for 
microstructural  characterization  and  mechanical  testing  under  uniaxial  and  stress/strain  controlled  cyclic 
loading.  The  crystal  plasticity  model  for  bcc  materials  uses  a  thermally  activated  energy  theory  for  plastic 
flow,  self  and  latent  hardening,  kinematic  hardening,  as  well  as  yield  point  phenomena.  Material 
parameters  are  calibrated  from  experiments  using  a  genetic  algorithm  based  minimization  process.  The 
computational  model  is  validated  with  experiments  on  stress  and  strain  controlled  cyclic  loading.  The 
effect  of  grain  orientation  distributions  and  overall  loading  conditions  on  the  evolution  of  microstructural 
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Figure  8:  Comparison  of  experimental  and  simulated  stress -strain  plots  for  (a)  strain  controlled  tests  and 
(h)  stress-controlled  test. 
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4. 


RELATION  TO  CURRENT  STATE  OF  KNOWLEDGE 


Within  the  framework  of  modeling  plastic  deformation,  two  approaches  viz.  the  continuum  and  the  crystal 
plasticity  approaches,  have  emerged.  Motivated  by  experimental  observations,  the  continuum  approach 
introduces  phenomenological  models  with  scalar  or  tensor  internal  variables  whose  growth  is  determined 
by  appropriate  evolution  laws.  Various  continuum  models  proposed,  include  the  anisotropic  quadratic 
yield  function  of  Hill  [42],  yield  surface  vertex  models  with  deviation  from  normality  [43,44],  anisotropic 
yield  criteria  with  isotropic  plasticity  equivalent  concept  [45],  and  isotropic  elastic-viscoplastic  models 
[46,47].  Though  popular  for  their  relative  simplicity,  they  do  not  properly  represent  underlying  physical 
complexities  due  to  the  lack  of  appropriate  length  scales.  Higher  order  gradients  have  also  been 
introduced  in  expressions  for  flow  stress  and  yield  surfaces  to  incorporate  length  scale  effects  in  these 
continuum  theories  [46,47].  However  they  still  lack  accurate  representation  of  those  phenomena,  which 
are  strongly  affected  by  evolution  and  interactions  on  the  microscale. 

In  the  crystallographic  approach,  the  mechanical  response  of  polycrystalline  aggregates  are  deduced  from 
the  behavior  of  constituent  crystal  grains  with  specific  assumptions  about  their  interaction.  For  example, 
the  Taylor  model  [49]  assumes  that  each  grain  is  subject  to  the  same  uniform  deformation  as  the  entire 
aggregate.  In  this  model,  compatibility  is  satisfied  and  equilibrium  holds  in  each  grain,  but  equilibrium  is 
usually  violated  between  grains.  Using  the  principle  of  maximum  work.  Bishop  and  Hill  [49]  calculated 
the  yield  stress  based  on  Taylor's  assumption.  The  Taylor-Bishop-Hill  theories  have  been  extensively 
used  to  predict  textures  and  stress-strain  response  in  fee  and  bee  polycrystalline  metals.  Though  the 
Taylor  model  and  its  derivatives  are  reasonably  effective  in  predicting  gross  features  of  deformation,  a 
number  of  studies  (e.g.  [50,51])  have  predicted  sharper  textures  than  experiments.  This  is  attributed  to 
nonuniform  deformations  due  to  grain  interactions  in  real  polyerystals  [52].  Relaxed  constraints  methods 
[53]  have  been  proposed  to  incorporate  intergranular  interactions,  with  relaxation  of  the  requirements  of 
Taylor  model.  Experiments  on  large  grain  sheet  materials  by  Skalli  et.al.  [54]  have  shown  that  full 
constraint  and  relaxed  constraint  models  each  work  well  only  for  some  grain  orientations.  Other  important 
contributions  in  this  area  include  the  work  of  Iwakuma  and  Ncmat-Nasser  [55],  who  have  used  Hill's  self 
consistent  averaging  methods  to  develop  finite  rotation-small  strain  constitutive  relations  for  rate 
independent  and  rate  dependent  crystals.  Most  of  these  models  have  however  overlooked  the  relation 
between  plastic  instabilities  and  the  details  of  micro  structural  description  and  intergranular  interactions. 

Various  computational  studies  have  modeled  anisotropy  and  its  evolution  in  large  deformation  processes 
with  the  crystallographic  approach  (e.g.  Havner  [56],  Mathur  et  al.  [57],  Kalidindi  et.  ah  [58]).  A 
majority  of  these  models  are  based  on  variants  of  the  Taylor  assumption,  thereby  trivializing  grain 
interactions  and  morphological  effects  like  grain  shape,  size  and  arrangements.  While  some  models  have 
utilized  continuous  orientation  distribution  functions  (ODFs)  [59],  a  majority  of  them  do  not  account  for 
microstructural  interactions  and  morphology,  in  predicting  the  onset  and  growth  of  plastic  localization. 
The  role  of  grain  shape  and  grain  boundaries  on  localization  has  been  studied  by  Lemonds  et,  al.  [60]  for 
bicrystals  of  identical  geometry,  but  of  different  crystallographic  orientations  and  hardening.  They 
concluded  that  grain  boundaries  cause  nonuniform  deformation  and  that  the  formation  of  shear  bands  is 
strongly  affected  by  initial  misorientation  across  the  boundary  and  the  difference  in  latent  hardening. 
Nonuniform  deformation  behavior  of  polycrystals  with  inter-granular  interaction  has  been  investigated  by 
Harren  and  Asaro  [50]  using  finite  element  models  with  very  large  number  of  elements  (element  size  ~ 
1/100  grain  size).  The  importance  of  grain  boundaries  and  triple  point  junctions  in  slip  activity  and 
deformation  nonuniformities  is  brought  out  in  their  research.  Similar  simulations  have  also  been 
conducted  by  Kalidindi  et.  al.  [51]  and  Beaudoin  et,  at.  [61].  Panchan adees waran  and  Doherty  [62]  have 
emphasized  the  need  for  incorporation  of  actual  metallurgical  structures ,  such  as  grain  size,  shape  and 
orientation  in  models,  due  to  otherwise  poor  qualitative  agreement  between  texture  predictions  and 
experimental  measurements.  This  has  also  been  emphasized  by  Becker  and  Richmond  [63]  in  their  finite 
element  simulations  of  TEM  micrographs  for  predicting  flow  localization  and  failure  in  aluminum  alloys. 
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Since  poly  crystal  line  deformation  is  predominantly  three-dimensional  in  nature,  it  is  essential  that  the 
microstructural  models  be  developed  with  detailed  3D  information.  This  requirement  poses  major 
challenges  to  currently  available  metallography  practices*  Techniques  based  on  ultrasonics  or  its  variants* 
such  as  acoustic  microscopy  or  laser  ultrasonics  [64*65]  rely  on  good  reflection  properties  and  have 
limited  application  in  metals.  While  X-ray  based  computed  tomography  [66*67]  methods  are  widely  used 
in  3D  solid  mode!  generation*  they  are  generally  deficient  in  achieving  the  resolution  desired  for  detailed 
study  of  polycrystalline  metals.  Recently*  synchrotron  based  CT  technology  has  been  developed  to  yield 
tomographic  images  with  considerably  high  resolution  [68].  However*  this  method  is  still  not 
commercially  available  and  in  generally  is  quite  expensive. 

Many  metals  and  their  alloys  exhibit  strong  grain  size  effects  due  to  strong  inhomogeneous  deformation 
fields  in  the  grain  interior  [69],  The  conventional  homogeneous  constitutive  models  of  [51,6!]  etc*  fail  to 
accurately  predict  the  scale  dependent  behavior.  The  nonuniform  deformation  occurs  at  length  scales, 
which  for  plastic  deformation  are  associated  with  the  dislocation  structures.  This  causes  the  plastic  flow  to 
be  size  dependent  with  an  increasing  hardening  for  decreasing  grain  sizes  (Hall  Petch  effect)  [70].  The 
nonuniform  plastic  deformations  give  rise  to  the  development  of  geometrically  necessary  dislocation 
(GND)  which  arc  required  to  preserve  crystallographic  lattice  compatibility  and  cause  short  and  long 
range  effects  on  the  constitutive  behavior.  The  GND5s  act  as  obstacles  to  the  motion  of  gliding 
dislocations  which  carry  the  plastic  deformation*  also  known  as  statistically  stored  dislocations  (SSD).  Of 
the  total  dislocation  density,  the  GND  density  is  defined  as  the  minimum  density  required  to 
accommodate  a  given  strain  gradient  and  SSD  density  is  the  rest  [71]*  Various  strain  gradient  models 
have  been  proposed  for  capturing  size  scale  effects.  Some  of  these  approaches  have  introduced  intrinsic 
material  length  scale  in  plasticity  [72*73,74,75,76],  while  others  have  related  the  underlying  physics  to  the 
generation  and  distribution  of  dislocations  [77]  or  to  plastic  and  elastic  incompatibilities  [78,79,80]. 
However  some  of  the  shortcomings  of  these  models  are  due  to  2D  idealizations,  considerations  of  global 
strain  grain  instead  of  those  at  the  grain  level,  and  global  incompatibilities  through  the  curl  of  the 
deformation  gradient  tensor  rather  than  slip  system  incompatibilities*  Additionally*  the  influences  of  the 
GND’s  in  these  models  are  superposed  on  the  hardening  or  slip  resistance,  thus  ignoring  their  distinct 
contributions  compared  to  SSD’s.  In  a  more  physically  based  modeling  effort,  Evers  [81]  has  developed  a 
model  where  the  GND  densities  are  directly  related  to  the  gradients  of  crystallographic  slip.  In  this  model, 
the  slip  resistance  is  caused  by  the  short  range  interaction  of  all  SSDs  and  GNDs.  Furthermore,  any 
heterogeneity  in  the  GND  density  field  also  causes  a  long  range  influence  on  crystallographic  slip  through 
back  stress  contribution,  since  at  a  length  scale  that  is  several  orders  of  magnitude  larger  than  each 
individual  dislocation  size*  the  net  effect  of  the  GND  densities  does  not  vanish.  This  physically 
motivated  strain  gradient  plasticity  model  strain  gradient  model  shows  promise  and  will  be  experimented 
with  in  the  proposed  work* 

The  recent  years  have  seen  significant  efforts  in  modeling  cyclic  plasticity  and  fatigue  with  considerations 
of  microstructural  stress-strain  evolution.  On  the  micro-scale,  the  state  of  stress  in  a  polycrystalline 
aggregate  during  cyclic  loading  depends  on  grain  orientation  in  relation  to  its  neighbors  as  well  as  with 
respect  to  the  loading  direction  and  constitutive  properties  of  the  material.  Finite  element  calculations 
have  shown  that  depending  on  the  loading  conditions,  significant  gradients  of  stress  can  arise  within  a 
single  grain.  The  variation  of  stress  over  grains  in  a  polycrystalline  aggregate,  at  any  given  location  in  the 
continuum  depends  on  the  single  crystal  properties  of  the  material  as  well  as  the  local  crystallographic 
texture.  Dawson  and  coworkeis  [82]  have  investigated  the  mechanical  behavior  of  aluminum  alloy  in 
cyclic  loading  using  crystal  plasticity  based  FEM  simulations  of  crystalline  aggregates.  Up  to  1000  cycles 
are  simulated  in  their  experiments.  While  initial  orientations  are  assigned  by  sampling  orientation 
distribution  functions  from  measured  crystallographic  orientations,  grain  morphology  is  not  explicitly 
accounted  for.  Neutron-diffraction  is  used  to  measure  residual  lattice  strains*  Their  simulations  show  that 
under  constant  amplitude  strain  cycles,  the  macroscopic  stress  decays  with  increasing  cycles  with 
evolving  lattice  strains.  Miller  and  Dawson  [83]  have  considered  a  feature  based  'digital  material’  with 
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statistical  description  of  morphology,  lattice  orientation,  intra-grain  lattice  misorientation,  and  material 
constitutive  model  in  an  attempt  to  deploy  accelerated  insertion  of  materials.  McDowell  and  coworkers 
have  incorporated  a  crystal  plasticity  model  with  kinematic  hardening  in  ABAQUS  to  model  cyclic 
plasticity  in  high  cycle  fatigue  in  Ti-6A1-4V  [84].  Random  distributions  of  grains  with  hep  secondary  a 
and  bee  [J  in  a  lamellar  structure  are  modeled  and  the  microstructure  is  subjected  to  cyclic  loading. 
Simulations  of  si  mple  components  have  also  been  done  to  consider  fretting  fatigue. 

A  number  of  fatigue  models  considering  the  effects  of  microstmctural  variations  have  been  considered  in 
the  recent  years.  Chan  [85]  has  formulated  fatigue  crack  growth  laws  in  the  power  law  regime  in  terms 
of  material  properties  like  yield  stress  and  fatigue  ductility  and  micros  true  turn  I  parameters  like  dislocation 
cell  size  and  dislocation  barrier  spacing.  He  attributed  the  lack  of  strong  microstmctural  influence  on  the 
fatigue  crack  growth  due  to  increasing  yield  stress  and  fatigue  ductility  with  decreasing  dislocation  barrier 
spacing,  A  probabilistic  treatment  of  this  model  has  been  considered  in  [86],  where  randomness  in  yield 
stress,  fatigue  ductility  and  dislocation  barrier  spacing  have  been  correlated  to  the  variations  in  fatigue 
crack  growth  rate.  The  recent  research  of  Chan  [87]  is  extending  the  above  models  for  describing 
dislocation  structure  evolution,  fatigue  crack  initiation  and  growth  and  linking  various  scales.  A  fatigue 
crack  initiation  model  is  constructed  in  terms  of  microstmctural  parameters  like  grain  size,  slip  band 
spacing,  inclusion  size  crack  size  and  notch  size.  Try  on  and  Dey  [88]  has  proposed  a  similar  methodology 
to  aid  in  the  design  and  prediction  of  reliability  of  metallic  components  based  on  probabilistic 
microstmctural  fatigue  models.  The  variation  of  fatigue  life  is  predicted  as  a  function  of  the  statistical 
variation  of  the  microstructure  of  the  material. 

The  proposed  research  has  sufficient  uniqueness  and  the  potential  to  advance  the  state  of  the  art  in  the 
modeling  of  fatigue  failure  in  real  polyeryslallinc  microstructures.  The  fatigue  crack  initiation  and  growth 
models  will  directly  evolve  from  finite  element  simulation  of  the  microstructure.  A  novel  contribution  of 
this  research  will  be  in  the  development  of  robust  computational  and  constitutive  models  for  prediction  of 
localization  and  subsequent  crack  evolution.  Adaptivity  of  computational  models  to  accommodate  the 
material  evolution  is  a  necessary  tool  that  has  not  been  applied  extensively  to  this  class  of  problems. 


5.  OUTLINE  OF  THE  OVERALL  APPROACH 

The  system  of  computational  models  and  codes  to  be  developed  in  the  proposed  research  will  be  an 
ensemble  of  tasks  that  are  elaborated  below.  It  is  expected,  that  on  completion,  the  project  will  be  able  to 
deliver  a  comprehensive  set  of  computational  tools  that  can  predict  the  effect  of  material  microstructure 
on  cyclic  deformation  and  fatigue  crack  evolution.  It  is  important  to  note  that  components  of  the  overall 
framework  will  be  modular,  in  that  they  can  be  interchanged  with  more  appropriate  models  as  they 
evolve.  This  PTs  ongoing  collaborations  with  Professors  Mills  and  Fraser  at  OSU  and  Drs.  Dirniduk, 
Parthasarathy  and  Uchic  at  AFRL/MLLM  will  be  critical  in  the  development  and  validation  of  the 
computational  models  through  results  of  mechanical  testing  and  quantitative  metallography.  Letters  from 
these  collaborators  are  attached  in  the  appendix.  Major  developmental  tasks  are  summarized  below. 


5,1  Microstmctural  Computer  Model  from  OIM  Data  and  Characterization 

Accurate  mapping  of  polycrystalline  micros truclures  from  metal lographic  observations  to  the  finite 
element  model  is  essential  for  reliable  prediction  of  deformation  and  damage  evolution  leading  to  failure. 
Mean  values  of  characteristic  features  like  phase  volume  fraction,  grain  size,  shape  and  orientation 
distributions  are  not  adequate  in  this  respect,  especially  with  respect  to  predictions  of  localization,  crack 
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initiation  and  growth.  A  combination  of  those  features  set  the  grain-level  kinematics  of  the  material  and 
controls  the  local  driving  forces  for  energy  partitioning,  defect  formation,  and  ultimately,  failure. 
Capturing  details  of  variabilities  of  key  3D  microstructural  features  is  important  to  enhance  the  predictive 
capability  of  the  computational  model.  This  requirement  poses  major  challenges  to  currently  available 
metallography  practices.  In  this  research,  the  methods  described  in  section  3.2  and  [2]  will  be  advanced  to 
create  equivalent  3-D  microstructural  images  from  metal lograp hie  analyses  data  as  described  next, 

A.  Sectioning-Bascd  3D  Mierostriicture:  flic  newly  developed  automated  approach  of  characterizing 
3D  microstructure  using  a  Dual  Beam  Focused  Ion  Beam  (FtB)-SEM  system  will  be  advanced  to  create  a 
seamless  transition  from  3-D  orientation  data  to  computer  models  of  model  grain  and  sub-grain  structures. 
The  microscope  is  capable  of  highly-localized  micro-machining  and  ion  imaging  using  the  FIB  column, 
and  non-destructive  high -resolution  imaging  and  Electron  Back-Scatter  Diffraction  (EBSD)  using  the 
SEM  column.  Serial  sectioning  of  the  specimen  by  the  FIB  is  followed  by  the  EBSD  system  to  determine 
grain  section  orientations.  The  Micro- Mcsher  program  has  been  developed  to  post-process  the  OIM  data 
and  generate  grain  boundary  structures  from  differences  in  grain  orientation  data  [2],  It  approximates 
spatially  complex  boundaries  with  line  segments  with  a  user  prescribed  error-per-unit- length  algorithm. 
This  program  will  be  refined  and  tested  in  this  research  for  a  variety  of  different  microstructures.  A 
sequence  of  the  FIB/EBSD/Micro-Mesher  generated  sections  (see  figure  3b)  will  be  stacked  using  3D 
visualization  software  like  Slicer-Dieer  [89]  to  construct  the  3D  microstructures.  The  2D  polycrystalline 
sections  and  the  3D  computer  generated  micro-regions  will  be  processed  for  microstructural 
characterization.  Important  microstructural  parameters  like  distributions  of  grain  size  and  shape, 
orientation  distribution  functions  (ODF)  and  other  stereological  characteristics  will  be  compared  for  2D 
and  3D  microstructures. 

B,  Statistically  Equivalent  Computed  Microstructure:  In  an  alternate  approach,  an  attempt  will  be 
made  to  generate  statistically  equivalent  3D  microstructures  from  surface  OIM  data  of  a  representative 
micro-region.  The  proposed  method  will  design  an  algorithm  to  generate  an  optimal  distribution  of  points 
for  3D  tessellation  of  the  micro-region,  based  on  these  points  or  seeds  (figure  9).  The  convergence  of  the 
optimal  point  distribution  algorithm  will  be  established  by  comparing  the  distributions  of  microstructural 
parameters  on  the  surface  of  the  tessellated  region  with  that  from  the  surface  OIM  data.  An  iterative 
algorithm  will  be  pursued,  till  the  convergence  criteria  are  satisfactorily  met  for  important  distribution 
functions.  This  method  is  expected  to  be  combined  with  established  stereological  methods  [90-94]  which 
estimate  the  number  density  and  size  of  3D  objects  from  density  and  size  of  2D  sections.  The  PI  has 
compared  predictions  of  some  stereological  methods  with  serial  sectioning  methods  to  examine  their 
effectiveness  in  [23,24].  This  experience  will  be  put  to  use  in  the  above  method.  The  experimental 
method  in  task  A  will  be  used  to  validate  the  proposed  method  in  task  B. 


Figure  9:  A  point  based  3D  tessellation  for  generating  a  statistically  equivalent  micro-region 
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5.2  Constitutive  Model  Development  and  Validation 

This  task  will  focus  on  the  incorporation,  development  and  validation  of  crystal  plasticity  constitutive 
models  for  specific  applications  of  polycrystalline  aggregates,  e.g.  Nickel  base  superalloys  like  NijAl 
(phases  with  fee  structure)  and  Ti  alloys  Ti-6  A1  (hep)  and  Ti  6242  (colonies  with  hep  with  bcc  phase 
structure).  Both  classes  of  materials  may  contain  more  than  one  phase  and  are  similar  in  size  scale.  The 
Ni  grains  are  essentially  equiaxed  with  particles  inside.  The  Ti  grains  may  have  colonies  of  similarly 
oriented  laths  of  a  (hep)  and  P  (bcc),  with  individual  laths  having  different  crystallographic  orientations. 
The  y/y*  boundary  in  Ni  is  relatively  low  energy  and  can  be  coherent.  In  Ti  the  a/p  interface  has  an 
orientation  relation,  but  is  not  as  low  energy  as  the  y/y’  boundary.  Important  aspects  of  this  work  are 
discussed  below. 


A.  Framework  of  the  Constitutive  Model:  The  primary  aim  of  this  task  will  be  to  examine  the 
effectiveness  of  constitutive  models  through  finite  element  modeling  of  single  crystal  and  polycrystalline 
behavior  using  crystal  plasticity  models,  implemented  in  the  finite  element  code  ABAQUS  through  the 
user  interface  UMAT.  Similar  implementation  of  experimentally  validated  models  has  been  accomplished 
by  the  Pi’s  group  in  [1,40,41].  Rate  dependent  crystal  plasticity  models  [95,96,9,98,99,100]  relate  the 

-  (X 

plastic  shear  strain  rate  Y  of  each  slip  systems  of  the  crystal  lattice  to  the  effective  shear  stress  r"  and 
the  slip  system  deformation  resistance  sa  through  a  visco-plastic  power  law  relation  as. 


where  m  characterizes  the  material  rate-sensitivity  and 


is  the  resolved  effective 


stress,  which  is  the  driving  force  for  the  dislocation  motion  on  slip  system  a,  is  the  resolved  back 

stress;  t  is  the  Schmid  stress,  &a  is  the  part  of  the  critical  slip  resistance  due  to  the  athermal  obstacles 

to  slip.  A(7  is  the  activation  free  energy  required  to  overcome  the  obstacles  to  slip  without  (he  aid  of  an 
applied  shear  stress;  k  is  the  Boltzmann’s  constant  and  T  is  the  temperature,  The  slip  system  resistance  is 

assumed  to  evolve  as:  sa  =  ,  where  huii  describes  the  rate  of  strain  hardening  on  slip  system 

a  due  to  shearing  on  the  slip  system  /?.  Kinematic  hardening  is  essential  for  modeling  cyclic 
deformation,  for  which,  an  evolution  law  for  the  slip  system  back  stress  has  been  proposed  by  Armstrong 
and  Frederick  [101]  as:  t"jn  =  (  ■  n"]hi“  ,  where  ^  <8>  n®  +  11"  ®m*).  The 

evolution  of  kinematic  hardening  completely  characterizes  reversible  non-proportional  cyclic  stress-strain 
behavior. 


B.  Microstruciural  Finite  Element  Models  of  Single  and  Polvcrvstals,  Colonies  etc.:  In  this  task, 
crystal  plasticity  parameters  will  be  calibrated  from  experiments  performed  by  Mills,  Dimiduk  (see 
letters)  by  inverse  methods.  Explicit  models  of  microstructural  regions  with  constituent  phases  will  be 
modeled  and  loaded  for  specific  slip  activation  in  accordance  with  the  experimental  studies.  For  example, 
in  a  Ti-6242  colony  specimen,  lath  regions  of  a  (hep)  and  p  (bcc)  phases  with  the  Burger’s  orientation 
relation  will  be  modeled  in  FEM  to  calibrate  slip  system  parameters  for  the  individual  phases  (see  figure 
10).  It  is  computationally  exhaustive  to  use  these  explicit  two-phase  models  to  represent  grains/colonies 
polycrystalline  aggregate  modeling.  Consequently  homogenized  effective  models  with  reduced  degrees  of 
freedom  will  be  developed.  In  these  models,  the  effect  of  the  explicit  phases  can  be  represented  by 


18 


corresponding  slip  systems  and  other  averaged  parameters  like  volume  fraction.  The  effective  models  will 
be  tested  for  polycrystalline  aggregates,  constructed  from  equivalent  orientation  distributions  of  those 
measured  via  el ectron-backsca tiered  diffraction  (EBSD)  techniques.  Results  of  the  polycrystaHine  model 
will  then  be  compared  with  experimentally  measured  constant  strain  rate  and  creep  tests  for  different 
orientations 
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Figure  10:  (a)  A  computational  model  of  transformed  fi  colony  showing  a  and  J3  phases  corresponding  to 
a3  basal  loading ,  (b)  a  comparison  of  experiments  and  simulated  results, 

C  Tension-Compression  Asymmetry:  Work  by  Mills  has  demonstrated  that  creep  rates  in 
polycrystalline  Ti-6242  are  significantly  larger  in  tension  than  in  compression.  A  question  to  be  answered 
in  this  research  is  whether  the  asymmetry  is  due  to  the  presence  of  (3  in  the  two  phase  alloys,  or  whether 
the  asymmetry  is  a  feature  of  the  a  phase  of  Ti  alloys  (perhaps  due  to  a  non-planar  dislocation  core 
configuration*  To  this  end.  Mills  will  be  conducting  tension  and  compression  testing  of  single  crystals  of 
Ti-6AL  This  information  will  then  be  used  by  to  develop  a  more  complete  description  of  the  deformation 
of  the  a  phase  which  may  include  this  asymmetry, 

IK  Grain  Size  Effects  in  Crystal  Plasticity  [Models  :  Many  metals  and  their  alloys  e.g.  Ti  6242,  exhibit 
strong  grain  size  effects  due  to  strong  inhomogeneous  deformation  fields  in  the  grain  interior  [69],  There 
have  been  few  direct  measurements  of  grain  size  effects  in  Ti  alloys,  and  no  published  studies  of  the 
effect  of  grain  size  on  creep  response.  Mills  at  OSU  is  pursuing  such  studies  on  single  phase  Ti-6A1, 
which  has  the  advantage  of  offering  a  relatively  “clean”  measurement  of  grain  size  effects,  unlike  studies 
of  two-phase  a-p  alloys  where  the  shape  and  dimension  of  the  individual  phases  can  obscure  the  grain 
size  effects.  These  results  will  be  used  to  calibrate  the  crystal  plasticity  model,  thereby  including  scale 
effects  directly.  Due  to  the  planar  nature  of  slip  in  these  short-range-ordered  alloys,  extended  pile-ups, 
local  stress  concentrations  and  grain -to -grain  “communication”  of  slip  is  observed  experi  mentally  using 
SEM  and  TEM, 

Consistent  with  the  discussion  in  section  4,  this  work  will  explore  size  dependent  models  in  which  strain 
gradient  related  production  of  GND’s  enable  the  prediction  of  size  dependent  behavior  for  intragranular 
and  global  strain  gradients  [81].  In  this  model  the  evolution  of  SSD’s  during  crystallographic  slip 
increases  the  number  of  short  range  interactions,  which  leads  to  isotropic  hardening.  The  SSD  density  is 
affected  by  generation  and  annihilation  of  SSD’s  on  each  slip  system  and  is  expressed  in  the  form: 
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due  to  a  balance  of  accumulation  and  annihilation  of  dislocation 

density  rate.  The  geometrically  necessary  edge  and  screw  dislocations  can  be  obtained  from  relations  to 
the  gradients  of  plastic  slip  corresponding  to  the  undeformed  configuration  as  [76,102]: 

+— •/>£)  .  Any  gradient  normal  to  the  glide  plane  do 

.  f  .  .  c 

* - v - J  -  V  - '  •  J 

edge  screw 

not  affect  the  lattice,  A  nonuni  form  GND  distribution  results  in  additional  residual  stresses  and 
consequently  a  back  stress  term  is  entered  in  the  crystal  plasticity  formulation.  This  counteracts  the  local 
resolved  shear  stress  and  obstructs  crystallographic  slip  through  long  range  kinematic  hardening.  Such 
models  will  be  incorporated  and  will  be  tested  with  experimental  results.  Parameters  in  the  size  dependent 
hardening  model  will  be  calibrated  from  experiments  on  polycrystals  of  varying  sizes. 


5.3  Fatigue  Crack  Evolution  with  Cohesive  Zone  Models 


Fatigue  crack  evolution  in  the  microstructure  takes  place  in  multiple  stages.  In  the  first  phase,  plastic  flow 
concentration  due  to  crystallographic  slip  and  interaction  between  grains  lead  to  intense  micro-shear  band 
formation  in  the  grains.  The  second  is  the  crack  nucleation  phase,  in  which  damage  accumulates  to  form  a 
crack,  while  the  third  stage  is  the  small  crack  growth  phase,  where  the  crack  size  is  on  the  order  of  the 
microstructure.  The  final  stage  is  a  long  crack  growth  phase,  where  the  crack  size  is  large  compared  to  the 
microstructure.  Each  stage  is  driven  by  different  mechanisms  and  will  be  distinctly  modeled. 


The  fatigue  mechanisms  will  be  explicitly  modeled  in  this  task  for  polycry stalline  microstructures  to 
examine  the  effects  of  local  morphology.  Fatigue  crack  growth  is  the  result  of  the  competition  of  material 
separation  processes  near  the  crack  tip  and  the  deformation  processes  in  the  surrounding  material.  An 
irreversible  cohesive  zone  traction-separation  law,  which  represents  the  physical  processes  of  material 
deterioration,  will  be  used  to  model  crack  initiations  and  growth  in  the  material  microstructure.  A  variety 
of  cohesive  zone  models  have  been  proposed  in  literature  [103,104,105].  In  this  work,  an  irreversible 
CZM  of  the  form  developed  in  by  [106]  will  be  used  to  model  fatigue  crack  growth. 


Tn 

Tt 
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This  model  incorporates  loading  and  unloading  conditions,  as  well  as  effects  of  accumulation  of  damage 
during  cyclic  loading  on  cohesive  zone  model  parameters.  The  traction-separation  behavior  is  irreversible 
and  history  dependent,  due  to  the  fact  that  tractions  depend  on  the  current  state  of  damage,  and  hence  on 
the  loading  history.  Parameters  in  the  cohesive  zone  models  will  be  calibrated  using  inverse  methods 
using  observations  from  fatigue  experiments,  consisting  of  mechanical  testing  and  OIM  scans  of  selected 
regions.  The  mode!  applies  to  material  systems  where  crack  bridging  dominates  the  crack  growth 
behavior,  and  no  energy  dissipation  takes  place  outside  of  this  zone.  The  cohesive  zone  model  will  be 
tested  extensively  for  effectiveness  with  intra-,  inter-  or  grain  boundary  failure.  Figure  1 1  shows  the 
results  of  a  preliminary  polycrystalline  model  with  the  cohesive  zone  interface. 
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Figu  re  1 1 :  (a)  A  trad  ion-separation  law  in  the  cohesive  zone  model,  (h)  stress-strain  behavior  in  a  single 
loading-unloading  cycle. 


5.4  A  Novel  Grain-Based  Finite  Element  Model  (GBFEM)  for  Finite  Deformation 

General  purpose  finite  element  codes  like  ABAQUS  are  severely  limited  in  modeling  large 
microstructural  regions  of  polycrystalline  aggregates  due  to  prohibitively  high  computing  requirements. 
The  computing  demands  are  even  intensified  when  numerous  time  steps  are  required  for  modeling  cyclic 
deformation.  It  is  inevitable  that  specialized  models  of  computational  mechanics  endowed  with  adaptive 
capabilities  can  bring  considerable  payoffs  in  this  regard.  In  the  proposed  task,  a  special  grain  based 
adaptive  finite  element  model  (GBFEM)  will  be  developed  for  accurate  and  efficient  analyses  and 
overcome  limitations  of  conventional  FE  models  in  fatigue  life  prediction.  The  model  will  incorporate  the 
real  material  morphology  from  task  5.1.  Rate  dependent  crystal  plasticity  constitutive  models  and 
microstructural  cohesive  zone  damage,  as  discussed  in  tasks  5.2  and  5.3,  will  be  included  in  this  large 
deformation  GBFEM.  This  development  will  follow  the  highly  successful  Voronoi  Cell  FEM  (VCFEM) 
[25-31]  that  has  been  proven  to  offer  a  tremendous  advantage  in  accurate  and  efficient,  image-based 
multi-scale  modeling  of  nonuniform  composite  and  porous  materials.  Early  formulation  and  modeling  of 
GBFEM  have  shown  significant  promise  in  modeling  polycrystalline  materials  (see  figure  12). 


Figure  12:  A  preliminary  attempt  with  GBFEM. 

The  GBFEM  formulation  will  be  based  on  the  enhanced  strain  formulations  that  have  been  successful  in 
modeling  non-homogeneous  and  localized  strain  fields.  Each  element  in  GBFEM  will  manifest  a  single 
grain  or  crystal.  Plastic  flow  concentration  due  to  crystallographic  slip  and  interaction  between  grains  is  a 
source  of  micro-shear  band  formation.  GBFEM  elements  will  be  selectively  enriched  adaptively  to 
differentiate  between  regions  of  high  and  low  deformation  gradients  and  impart  preferential  resolutions  to 
sub-domains.  Adaptive  methodologies  will  be  instrumental  in  identifying  regions  of  evolving  localized 
deformation  and  preferential  resolutions  will  be  imparted  by  local  enrichment  of  sub-grain  domains.  In 
addition,  intra-,  inter-  and  grain  boundary  failure  will  be  modeled  by  using  the  cohesive  zone  models  in 
GBFEM  as  well.  Incremental  adaptivity  will  be  necessary  for  modeling  the  evolution  of  fracture  and 
bridging  nascent  surfaces. 
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A  necessary  ingredient  for  the  success  of  this  computationally  intensive  research  is  the  incorporation  of 
advanced  high  performance  computing  techniques.  Effective  use  of  high  performance  computing  will  be 
explored  to  reduce  the  overall  execution  time.  Novel  computational  frameworks  and  programming 
environments  will  be  developed  to  support  efficient  multi-scale  simulations  on  different  high  performance 
computing  platforms.  Parallel  programs  will  conform  to  the  latest  message-passing  interface  (MPI  and 
OpenMP)  standard  for  portability,  implementation  focus  will  be  on  aspects  of  (a)  Parallel  scalability,  (b) 
Mapping  and  scheduling,  (c)  (c)  Run-time  environment:  (d)  Parallel  I/O: 


5*5  Multi-Time  Scaling  Algorithm  for  Modeling  Long  Range  Cyclic  Loading 

Modeling  cyclic  deformation  using  conventional  time  increments  can  be  an  exorbitant  task  for  crystal 
plasticity  computations.  Most  simulations  performed  with  3D  crystal  plasticity  are  in  the  range  of  100 
cycles  [82,83,84]  and  the  results  are  subsequently  extrapolated  for  fatigue  predictions.  In  modeling 
fatigue,  it  is  however  desirable  to  conduct  simulations  for  a  significantly  high  number  of  cycles  to  reach 
local  states  of  damage  initiation  and  growth.  Conventional  methods  of  lime  integration  using  semi- 
discretization  present  numerous  challenges  due  to  the  variation  in  time  scales  ranging  from  the  scale  of 
the  entire  process  to  the  time  resolution  required  by  the  damage  evolution  and  crack  propagation 

For  this  purpose,  a  multiple  time  scaling  methodology  will  be  developed  in  this  task,  for  compression  and 
rarefaction  of  time  scales.  This  will  enable  large  time  scale  homogenization  in  relatively  benign  periods  of 
deformation,  to  very  fine  time  scale  simulations  in  temporal  regions  of  evolving  localization  or  damage. 
The  multi-time  scaling  will  based  on  homogenization  with  the  asymptotic  expansion  method  in  the  time 
domain  for  poly  crystals.  This  basic  idea  of  multi-time  scaling  has  been  developed  in  Asai  and  Ghosh  [4], 
where  a  fast  time  scale  t,  and  a  slow  time  scale  t  are  related  by  the  relation  *  =  ^ ,  where  £  is  a  small 
number.  Variables  in  the  crystal  plasticity  laws  can  be  written  using  the  double  temporal  scale  asymptotic 
expansion  is  of  the  form  u£  =w°(x:,fjr)+ftfl{jtTrTr)+^(£2)  -  The  various  order  terms  can  then  be 
uncoupled  and  averaged  to  transcend  time  scales.  Figure  13  shows  a  comparison  of  a  cyclic  deformation 
problem  solved  by  (a)  a  detailed  single  time  scale  time  integration  scheme  called  reference  solution,  fb)  a 
homogenized  solution  using  a  compressed  slow  scale  and  (c)  a  region  of  de- compress  ion  at  a  critical 
region  to  revert  to  the  detailed  time  resolution.  The  comparison  with  the  reference  solution  shows 
excellent  accuracy  while  the  efficiency  gained  through  time  compression  can  be  enormous.  This 
promising  technique  will  be  developed  to  complete  fruition  in  this  task. 


Figure  13:  Sir  ess -time  plot  in  a  cyclic  deformation  simulation  sho  wing  different  methods  of  time 
integration. 
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In  summary,  the  tasks  proposed  for  this  complex  set  of  problems  are  highly  challenging.  The  PI  believes 
that  his  prior  and  current  accomplishments  in  related  fields,  and  the  collaboration  with  experts  in 
complimentary  Helds  will  give  him  a  strong  impetus  to  satisfactorily  execute  this  ambitious  project. 


6.  IMPACT  ON  AIR  FORCE  SYSTEMS 

Assessment  and  enhancement  of  life  of  structural  materials  in  aerospace  structural  components, 
especially  in  aircraft  engines  is  currently  a  priority  at  the  US  Air  Force  and  other  agencies.  A 
number  of  multi-year  major  funding  opportunities  like  the  DARPA  AIM  and  Prognosis 
initiatives  have  been  announced  in  the  recent  years  to  address  this  very  important  but  challenging 
problem.  The  scope  of  the  proposed  research  project  fits  well  into  the  grand  scheme  of  objectives. 
Computational  mechanics  has  traditionally  played  and  will  continue  to  play  an  important  role  in  the 
development  of  novel  materials  and  systems.  The  proposed  project  aims  to  advance  the  state  of  the  art  in 
computational  mechanics  for  the  specific  objective  of  better  prediction  of  matcrial/structural  life  and 
reliability.  It  is  expected  that  the  approaches  for  reliable  analysis  and  design  will  be  established  and 
technology  will  be  transferred  to  the  engine  design  community.  The  model  and  technology  will  benefit 
various  applications  of  metallic  materials  with  high  structural  integrity  for  aerospace  systems, 

7.  DESCRIPTION  OF  ORGANIZATION 

Information  on  the  Department  of  Mechanical  Engineering  at  the  Ohio  State  University  is 
available  on  the  web  [103].  OSU  has  had  extensive  contracting  experience  with  the  government. 

8.  OTHER  AGENCIES  RECEIVING  THIS  PROPOSAL 

No  other  funding  agency  will  receive  this  proposal 

9.  IMPACT  ON  ENVIRONMENT 

The  proposed  effort  will  have  no  impact  on  water,  atmosphere,  natural  resources  and  human 
resources. 

10.  PROPOSED  PROJECT  SCHEDULE 
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JUSTIFICATION  OF  BUDGET 


•  S.  Ghosh 

Support  for  1  summer  month  is  being  requested.  The  PI  will  coordinate  the  entire  collaborative  effort  and 
be  involved  in  the  development  of  computational  models,  in  the  planning  and  execution  of  experiments, 
and  be  the  link  between  all  participants  in  this  research.  He  will  be  responsible  for  all  deliverables 
resulting  from  this  research,  including  papers,  and  project  reports. 

•  Graduate  Students 

Support  for  2  Ph.D.  students  has  been  requested.  One  student  will  be  responsible  for  major  developments 
of  the  computational  models  and  algorithms  and  material  characterization.  The  other  student  will  be 
mainly  responsible  for  the  validation  with  experimental  results  and  probabilistic  modeling.  The  students 
will  also  help  with  the  experiments  at  AFRL  and  OSU-MSE. 

•  Travel 

Monthly  meetings  with  the  collaborators  at  AFRL  are  planned.  In  addition,  students  will  also  travel  to 
AFRL  to  conduct  experiments,  and  stay  there  for  extended  periods  as  required.  Two  national  meeting  will 
be  attended  every  year  to  present  the  developments. 

•  Services  at  No  Cost  to  AFOSR 

Dr.  Dennis  Dimiduk  and  his  research  group  will  participate  in  this  research  at  no  additional  cost  to 
AFOSR.  Specifically  a  student  will  be  supported  under  the  STW-21  initiative  to  provide  experimental 
assistance  to  the  program  under  their  supervision.  In  addition,  Professor  M.J.  Mills,  who  is  supported  by 
AFOSR  will  work  with  the  PI  for  experimental  validation  and  constitutive  modeling. 


12  DESCRIPTION  OF  AVAILABLE  FACILITIES 

Computer  code  developments  will  be  performed  on  workstations  in  the  Pi's  Computational  Research 
Laboratory  at  the  Ohio  State  University.  A  high-end  computer  workstation  is  being  requested  with  this 
proposal  for  code  development  and  debugging.  Large-scale  simulations  will  also  be  done  on  the  SGI 
Origin,  CRAY  SV1  and  Beowulf  cluster  at  the  Ohio  supercomputer  center. 

Experimental  component  of  the  proposed  research  will  be  conducted  in  the  Materials  Testing  Air  Force 
Research  Laboratory  Materials  and  Manufacturing  Directorate.  Facilities  at  AFRL  will  be  available  for 
the  program. 

12  LETTERS  OF  INTENT 

The  letters  of  intent  from  collaborators  are  included  near  the  end  of  this  proposal  . 

1.  Dr.  D.  Dimiduk ,  Materials  and  Manufacturing  Directorate,  AF  Research  Laboratory 

2.  Professor  M.  J.  Mills,  Materials  Science  &  Engineering,  The  Ohio  State  University 

3.  Professor  Hamish  Fraser,  Materials  Science  &  Engineering,  The  Ohio  State  University 


24 


REFERENCES 


1.  V.  Hasija,  S.  Ghosh,  MJ.  Mills  and  D.S.  Joseph,  '"Modeling  Deformation  and  Creep  in  T1-6AI 
alloys  with  Experimental  Validation”,  Acta  Materially  VoL  51,  pp.  4533-4549,  2003. 

2.  M.  Groeber,  13.  Haley,  M.  Uehic  and  S.  Ghosh,  Microslructural  Characterization  using  3-D 
Orientation  Data  Collected  by  an  Automated  FIB-EBSD  system.  Proceedings  of  NUM1FORM  2004 ,  A  IP 
Press,  S.  Ghosh,  J.  Castro  and  J.K.  Lee  (eds.),  June  2004. 

3.  ABAOUS.  http:/ wwu.hks.com. 

4.  M.  Asai  and  S.  Ghosh,  Multiple  time  scale  modeling  for  cyclic  deformation  with  crystal 
plasticity.  Proceedings  of  NU Ml FORM  2004 ,  A  IP  Press,  S.  Ghosh,  J.  Castro  and  J.K.  Lee  (eds.),  June 
2004. 

5.  D.B.  Miracle  and  B.  Maruyama,  Metal  matrix  composites  for  space  systems:  current  uses  and 
future  opportunities,  Proc.  Of  IHPRPT  Meeting,  San  Diego,  Apr  2000. 

6.  D.B.  Miracle  and  B.  Maruyama,  Metal  matrix  composites  for  space  systems:  current  uses  and 
future  opportunities,  Proc .  Of  IHPRPT  Meeting,  San  Diego,  Apr  2000. 

7.  L.F.  Coffin,  Fatigue,  Ann.  Rev.  Mails.  Sci. ,  2,  (1973),  3 13-. 

8.  C.  M,  Laird,  Cyclic  deformation  of  metals  and  alloys.  Treat  Mails.  Sci.  Tech.,  R.J.  Arsenault 
(ed),  6,(1975),  10I-. 

9.  N.A.  Fleck,  KJ.  Kang,  and  M.F,  Ashby,  The  cyclic  properties  of  engineering  materials,  Acta 
Metal  et  Mat  ,  42,  (1994),  365-. 

10.  T.M.  Hashimoto,  M.S.  Pereira,  Fatigue  life  studies  in  carbon  dual-phase  steels,  Int  J.  Fatigue , 
Vol.18,  No.B,  (1996),  529-. 

11.  S.  Suresh,  Fatigue  of  Materials,  Cambridge  University  Press ,  Cambridge  ,  England,  (1991). 

12.  P.C.  Paris,  The  fracture  mechanics  approach  to  fatigue,  Fatigue- An  Interdisciplinary  Approach, 
Syracuse  University  Press,  (1964). 

13.  D,C.  Lord  and  L.F.  Coffin,  Low  cycle  fatigue  hold  time  behavior  of  cast  Rene  80,  Metal.  Tram., 
4,(1973),  1 647-. 

14.  H.  Tsuji  and  T.  Kondo,  Strain-time  effects  in  low  cycle  fatigue  of  Nickel -based  heat  resistant 
alloys  at  high  temperature.  Jour .  Nucl .  Mater.,  190,  (1987),  259-. 

15.  S.D.  Antolovich,  S.  Liu  and  R.  Baur,  Low  cycle  fatigue  behavior  of  Rene  60  at  elevated 
temperature,  Metal.  Trans.,  12  A,  (1981),  473-. 

16.  J.S.  Kallend  and  G.J.  Davies,  The  elastic  and  plastic  anisotropy  of  cold-rolled  sheets  of  copper, 
gliding  metal  and  a  brass,  Jour.  Inst.  Metals ,  Voh  5,  pp  257-260,  I97L 


25 


17,  D.  Juul  Jensen  and  N.  Hansen,  Relations  between  texture  and  flow  stress  in  commercially  pure 
aluminum,  in  Constitutive  Relations  and  their  Physical  Basis ,  8th  Riso  International  Symposium  on 
Metallurgy  and  Materials  Science,  S,l.  Andersen  et.  al.  (eds.),  pp  353-360,  1987. 

18,  A.  Korbel  and  P,  Martin,  Microscopic  versus  macroscopic  aspect  of  shear  band  deformation,  Acta 
Metal!,  et  Mater.,  Vol.  38,  pp  1369-1380,  1990. 

19,  R.Q.  Chu,  Z.  Cai,  S.X.  Li,  Z.G.  Wang,  Fatigue  crack  initiation  and  propagation  in  an  ex-iron 
polycrystals.  Mater.  Sci.  &  Engng.,  A3 13  (2001)  64-68. 

20,  A.D.  Boyd-Lee.  Fatigue  crack  growth  resistant  microstructures  in  polycrystalline  Ni-base 
superalloys  for  aeroengines,  Inter.  Jour.  Fatigue,  21,  (1999)  393-. 

21,  V,  Brien,  B.  Decamps.  Low  cycle  fatigue  of  a  nickel  based  superalloy  at  high  temperature: 
deformation  microstructures,  Mater.  Sci.  &  Engng.,  A3 16  (2001)  18-. 


22.  M,  Li,  S.  Ghosh,  T.N.  Rouns,  H.  Weiland,  O.  Richmond  and  W.  Hunt,  Serial  sectioning  method 
in  the  construction  of  3-D  microslructures  for  particle  reinforced  MMCs,  Materials  Characterization, 
Vol.  41,  No.  2-3,  pp.  81-95,  1998. 

23.  M.  Li,  S.  Ghosh,  O.  Richmond,  H.  Weiland  and  T.N.  Rouns,  Three  dimensional  characterization 
and  modeling  of  particle  reinforced  MMC’s,  Part  I:  Quantitative  description  of  morphology,  Mater.  Sci. 
Engng.  A,  Mol  A265,  pp.  1 53- 1 73,  1 999. 

24.  M,  Li,  S.  Ghosh,  O.  Richmond,  H.  Weiland  and  T.N.  Rouns,  Three  dimensional  characterization 
and  modeling  of  particle  reinforced  MMC’s,  Part  II:  Damage  characterization.  Mater.  Sci.  Engng,  Vol. 
A266,  pp.221-240,  1999. 

25.  S.  Moorthy  and  S.  Ghosh,  Adaptivity  and  convergence  in  the  Voronoi  cell  finite  element  model 
for  analyzing  heterogeneous  materials.  Comp.  Meth  Appl.  Mech.  Engng.,  Vol.  185,  pp.  37-74,  2000. 

26.  S.  Moorthy  and  S.  Ghosh,  A  Model  for  Analysis  of  Arbitrary  Composite  and  Porous 
Microstructures  with  Voronoi  Cell  Finite  Elements,  Int.  Jour.  Numer,  Meth.  Engng.,  Vol,  39,  pp.  2363- 
2398,  1996. 

27.  S.  Ghosh  and  Y.  Liu,  Voronoi  cell  finite  clement  model  based  on  micropolar  theory  of  thermo¬ 
elasticity  for  heterogeneous  materials,  Int.  Jour.  Numer.  Meth.  Engng.,,  Vol.  38,  pp.  1361-1398,  1995. 

28.  S.  Ghosh  and  S.  Moorthy,  Particle  fracture  simulation  in  nonuniform  microstructures  of  metal- 
matrix  composites,  Act  a  Materialia,  Vol.  46,  No.  3,  pp.  965-982,  1998. 

29.  S.  Moorthy  and  S.  Ghosh,  A  Voronoi  cell  finite  element  model  for  particle  cracking  in  composite 
materials,  Comp.  Meth.  Appl.  Mech.  Engng.,  Vol.  15l,pp.  377-400,  1998. 

30.  S.  Ghosh,  Y.  Ling,  B.  S.  Majumdar  and  R.  Kim,  Interfacial  debonding  in  multiple  fiber- 
reinforced  composites,  Mech.  Materials,  Vol.  32  (10),  pp.  561  -591,  2000. 

31.  M.  Li,  S.  Ghosh  and  O.  Richmond,  An  experimental-computational  approach  to  the  investigation 
of  damage  evolution  in  discontinuously  reinforced  aluminum  matrix  composite,  Acta  Materialia,  Vol.  47 
No.  12,  pp.  3515-3532,  1999. 


26 


32.  S.  Ghosh,  K.  Lee  and  P.  Raghavan,  A  multi-level  computational  model  for  multiscale  damage 
analysis  in  composite  and  porous  materials.  Int.  Jour.  Solids  Struct..  Vol.  38  No.  14,  pp.  2335-2385, 2001. 

33.  P.  Raghavan,  S.  Moorthy,  S.  Ghosh  and  N.  Pagano,  Revisiting  the  composite  laminate  problem 
with  an  adaptive  multi-level  computational  model,  Comp.  Science  and  Technology,  (in  press). 

34.  K.  Lee,  S.  Moorthy  and  S.  Ghosh,  Multiple  scale  computational  model  for  damage  in  composite 
materials.  Comp.  Meth.  Appl.  Mech.  Engng.,  Vol.  172  No.  1-4,  pp.  175-201,  1999. 

35.  K.  Lee  and  S.  Ghosh,  'A  microstructure  based  numerical  method  for  constitutive  modeling  of 
composite  and  porous  materials.  Mater.  Sci.  Engng.  A,  Vol.  A272,  pp.  120-133,  1999. 

36.  S.  Ghosh,  K.  Lee  and  S.  Moorthy,  Two  Scale  Analysis  of  Heterogeneous  Elastic- Plastic  Materials 
with  Asymptotic  Homogenization  and  Voronoi  Cell  Finite  Element  Model,  Comp.  Meth.  Appl.  Mech. 
Engng. ,  Vol.  132,  pp.  63-116,  1996. 

37.  S.  Ghosh,  K.  Lee  and  S.  Moorthy,  Multiple  scale  analysis  of  heterogeneous  elastic  structures 
using  homogenization  theory  and  Voronoi  cell  finite  element  method,  Int.  Jour.  Solids  Struct.,  Vol.  32, 
No.  l,pp.  27-62,  1995. 

38.  P.  Raghavan,  S.  Li  and  S.  Ghosh,  "Two  Scale  Response  and  Damage  Modeling  of  Composite 
Materials”,  Finite  Elements  in  Analysis  and  Design,  (in  press). 

39.  P.  Raghavan  and  S.  Ghosh,  "  A  Continuum  Damage  Mechanics  Model  for  Unidirectional 
Composites  Undergoing  Interfacial  Debonding”  Mechanics  of  Materials,  (accepted for  publication). 

40.  C.L.  Xie,  M.  Grocber,  R.  Butler  and  S.  Ghosh,  Computational  Simulations  of  Tensile  and  Fatigue 
Tests  of  HSLA  Steel  Based  on  Experiments,  Proceedings  of  SAE  World  Congress,  2003. 

41 .  C.  Xie,  M.  Grocber  and  S.  Ghosh,  "  Modeling  Cyclic  Deformation  of  HSLA  Steels  Using  Crystal 
Plasticity  ",  ASME  Journal  of  Engineering.  Materials  and  Technology,  (accepted  for  publication). 

42.  R.  Hill,  The  Mathematical  Theory  of  Plasticity,  Clarendon  Press,  Oxford,  1950. 

43.  J.R.  Rice,  The  localization  of  plastic  deformation,  Theoretical  and  Applied  Mechanics,  W.T. 
Koiter  (cd.),  North-Holland,  pp  207-220,  1976. 

44.  A.  Needleman  and  J.R.  Rice,  Limits  to  ductility  set  by  plastic  flow  localization.  Mechanics  of 
Sheet  Metal  Forming,  D.P.  Koistinen  and  N.M.  Wang(eds-),  Plenum,  New  York,  pp  237-267,  1978. 

45.  A.P.  Karaftllis  and  M.C.  Boyce,  A  general  anisotropic  yield  criterion  using  bounds  and  a 
transformation  weighting  tensor,  Jour.  Mech.  Phys.  Solids,  Vol.  41,  pp  1859-1886,  1993. 

46.  S.B.  Brown,  K.H.  Kim  and  L.  Anand,  An  internal  variable  constitutive  model  for  hot  working  of 
metals,  Inter.  Jour.  Plasticity,  Vol.  5,  pp  95-130,  1989. 

47.  L.  Anand,  Elaslo- viscoplasticity:  Constitutive  modeling  and  deformation  processing.  Large 
Plastic  Deformations:  Fundamentals  of  Metal  Forming,  C.  Teodosiu,  J.L.  Raphanel,  F.  Sidoroff  (eds.), 
1993. 


27 


48.  G.I.  Taylor,  Plastic  strains  in  metals,  /  Inst.  Metals ,  Vol.  61,  pp  307-324,  1938. 

49.  J.F.W,  Bishop  and  R.  Hill,  The  theory  of  the  plastic  distortion  of  a  polycrystalline  aggregate 

under  combined  stresses,  Phil.  Mag.,  Vol.  42,  pp  414-427,  1951. 

50.  S.V.  Harren,  and  R.J.  Asaro,  Nonuni  form  deformations  in  polycrystals  and  aspects  of  the  validity 
of  the  Taylor  model, /.  Mech.  Phys.  Solids,  Vol.  37,  No.  2,  pp  191-232,  1989. 

51.  S.  R.  Kalidindi,  C.A.  Bronkhorst  and  L.  Anand,  On  the  accuracy  of  the  Taylor  assumption  in 
poly  crystal  line  plasticity,  Anisotropy  and  Localization  of  Plastic  Deformation ,  J.P.  Boehler  and  A.S. 
Khan  (eds.),  Elsevier  Science,  N Y,  pp  1 39- 1 42,  1991. 

52.  I.L.  Dillamore  and  H.  Katoh,  Jour.  Metal  Set,  Vol.  8,  pp  2] ,  1974. 

53.  G.R.  Canova,  U.F.  Kocks  and  J.J.  Jonas,  Theory  of  torsion  texture  development,  Acta  Metall.  et 

Mater.,  Vol.  37,  No.  2,  pp  191-232,  1984. 

54.  A.  Skalli,  R.  Fortunier,  R.  Fillit  and  J.H.  Driver,  Crystal  rotations  during  the  rolling  of  large 
grained  aluminum  sheet,  Acta  Metall.  et  Mater.,  Vol.  33,  pp  997-1007,  1985. 

55.  T.  Iwakuma  and  S.  Nemat-Nasser,  Finite  elastic-plastic  deformation  of  polyerystalline  materials, 
Proc.  R.  Soc.,  Vol.  A394,  pp  87,  1984. 

56.  K.S.  Havner,  Finite  Plastic  Deformation  of  Crystalline  Solids ,  Cambridge  University  Press,  NY, 
1992. 

57.  K.K.  Mathur,  P.R.  Dawson  and  U.F.  Kocks,  On  modeling  anisotropy  in  deformation  processes 
involving  textured  polycrystals  with  distorted  grain  shape,  Mech.  Mater.,  Vol.  10,  pp  183-202,  1990, 

58.  S.  R.  Kalidindi,  C.A.  Bronkhorst  and  L.  Anand,  Crystallographic  texture  evolution  in  bulk 
deformation  processing  of  fee  metals,  Jour.  Mech.  Phys.  Solids,  Vol.  40,  pp  537-569,  1992. 

59.  A,  Kumar  and  P.R.  Dawson,  The  simulation  of  texture  evolution  during  bulk  deformation 
processes  using  finite  elements  over  orientation  space.  Simulation  of  Materials  Processing;  Theory, 
Methods  and  Applications,  NUMIFORM'95,  S.F.  Shen  and  P.  Dawson  (eds.),  A.  A.  Balkema,  Rotterdam, 
pp  451-456,  1995. 

60.  J.  Lemonds,  R.J.  Asaro  and  A.  Needle  man,  A  numerical  study  of  localized  deformation  in 
bicrystals,  Mech.  Materials,  Vol.  4,  pp  363-373,  1985. 

61.  A.J.  Beaudoin,  H.  Mecking  and  U.F.  Kocks,  Development  of  local  shear  bands  and  orientation 
gradients  in  fee  polycrystals.  Simulation  of  Materials  Processing;  Theory,  Methods  and  Applications. 
NUMIFORM'95,  S.F.  Shen  and  P.  Dawson  (eds.),  A. A.  Balkema,  Rotterdam,  pp  225-230,  1995. 

62.  S.  Panchanadeeswaran  and  R.D.  Doherty,  Direct  observation  of  orientation  change  by  channel  die 
compression  in  polycrystalline  Aluminum-  Use  of  split  sample,  Scripta  Metallurgica  et  Mater.,  Vol.  28, 
pp  213-218,  1993. 

63.  R.  Becker  and  O.  Richmond,  Incorporation  of  microstructural  geometry  in  material  modeling, 
Modeling Simul.  Mater.  Sci.  Eng.,  Vol.  2,  No.  2,  pp  439-454,  1994. 


28 


64. 


S.K.  Burke,  S.  McK.  Cousland  and  C.M.  Scala,  Mater.  Forum  19,  Vol.  18,  pp.  85,  1994. 

65.  R.  G.  King,  and  P.M.  Delaney,  Mater.  Forum  19,  Vol.  18,  pp.  21,  1994. 

66.  M.  Berman,  L.M.  Bischof,  E.J.  Breen  and  G.M.  Pcden,  Mater.  Forum  19,  Vol.  18,  pp,  85,  1994. 

67.  J.Y.  Buffiere,  E.  Maire,  C.  Verdu,  et.  a!..  Mater.  Sci.  EngngA,  Vol.  633,  pp.  234,  1997. 

68.  R.N.  Yancey  and  G.Y.  Baaklini,  Proc.  AS  ME  Inter.  Gas  Turbine  Aeroengines  Congress, 
Cincinnati,  OH,  1993. 

69.  R.  Becker  and  S.  Panchanadeeswaran,  Effects  of  grain  interaction  on  deformation  and  texture, 
Acta  Metall.  Mater.,  Vol.  43,  pp.  2107-21 10,  1995. 

70.  A.  Needletnan  and  G.  Sevillano,  Preface  to  viewpoint  set  on:  geometrically  necessary 
dislocations  and  size  dependent  plasticity,  Scripta  Metallur.  Vol.  48,  109-1 1 1,  2003. 

71.  H.  Gao  and  Y.  Huang,  Scripta  Materialia,  Vol.  48,  pp.  1 13-1 18,  2003. 

72.  E.C.  Aifantis,  The  physics  of  plastic  deformations,  Inter.  Jour.  Plasticity ,  Vol.  3,  pp  21  1-247, 
1987. 

73.  H.M.  Zbib  and  E.C.  Aifantis,  On  the  localization  and  post-localization  behavior  of  plastic 
deformation,  1  and  II,  Res  Mechunica,  Vol.  23,  pp  261-291,  1988. 

74.  N.A.  Fleck  and  J.W.  Hutchinson,  Strain  gradient  plasticity,  Adv.  Appl.  Mech.,  Vol.  33,  pp  295- 
362,  1997. 

75.  J.  Y.  Shu,  N.A.  Fleck  and  W.E.  King,  Bicrystals  with  strain  gradient  effects,  Mat.  Res.  Soc. 
Symp.  Proc.,  Vol.  458,  pp  295-300,  1997. 

76.  Ashby,  M.F.,  The  deformation  of  plastically  non-homogeneous  materials.  Phil.  Mag.  21,  399- 
424,  1970. 

77.  H.  Gao,  Y.  Huang,  W.  Nix  and  J.W.  Hutchinson,  Mechanism  based  strain  gradient  plasticity-I, 
Theory.  Journal  of  Mech.  Phys.  Solids,  Vol.  47,  1239-1263,  1999. 

78.  H.  Dai,  Geometrically  necessary  dislocation  density  in  continuum  plasticity  theory;  FEM 
implementation  and  applications.  Ph.D.  thesis,  MIT,  1997. 

79.  Acharya  A.  and  Beaudoin  A.J,  2000.  Grain-size  effect  in  viscoplastic  polycrystals  at  moderate 
strains./  Mech.  Phys.  Solids  ,  Vol.  48,  2213-2230. 

80.  Arsenlis  A.  and  Parks  D.M.,  1999.  Crystallographic  aspects  of  geometrical ly-necessary  and 
statistically-stored  dislocation  density.  Acta  Materialia,  47:  (5),  1597-161 1. 

81.  L.  P.  Evers,  Strain  Gradient  crystal  plasticity  based  on  dislocation  densities,  Ph.D.  Dissertation, 
TU  Eindhoven,  Eindhoven,  The  Netherlands,  2003. 

82.  H.S,  Turkmen,  R.E.  Loge,  P.R.  Dawson  and  M.  Miller,  On  the  mechanical  behavior  of  AA  7075- 
T6  during  cyclic  loading,  Int.  Jour.  Fatigue,  Vol.  25,  pp.  267-281,  2003. 


29 


83.  M.  Miller  and  P.R.  Dawson,  Accelerating  the  evaluation  of  critical  mechanical  properties  of 
polyphase  alloys  by  integrated  simulation  and  experiment,  AIAA  Journal ,  (in  press). 

84.  Bennett  V.P.  and  McDowell  D.L.,  2003.  Polycrystal  orientation  distribution  effects  on  microslip 
in  high  cycle  fatigue.  InL  J.  Fatigue,  25, 27-39. 

85.  K.S.  Chan,  Scaling  laws  for  fatigue  crack  growth  of  large  cracks  in  steels.  Met.  Tram.  A,  Vol. 
24 A,  pp.  2473,  1993. 

86.  K.S.  Chan  and  T.Y.  Torng,  Probabilistic  treatment  of  microslruclura!  effects  on  fatigue  crack 
growth.  Jour.  Engng.  Mater.  Tech.,  Vol.  1 18, 379,  1996. 

87.  K.  S.  Chan,  Development  of  a  physically  based  methodology  for  predicting  material  variable  in 
fatigue  crack  initiation  and  growth  response,  AFOSR  Review,  2002. 

88.  R.G.  Tryon  and  A.  Dcy,  Reliability  based  computational  model  for  material  development  and 
structural  design.  Inter.  Jour.  Mater.  Product.  Tech.,  Vol  16.,  No.  4,  pp  333-357,  2001. 

89.  Pixotech  Corporation.  Slicer-Dicer,  Version  3.5,  User’s  Manual. 

90.  E.R.  Wei  bet,  Stereological  Methods:  Theoretical  Foundation,  Vol.  2,  Academic  Press,  London, 
1980. 

91.  J.C.  Russ,  Practical  Stereology,  Plenum  Press,  NY,  1986. 

92.  S.  A.  Saltykov,  Stereometric  Metallograph,  2nd  ed,  State  Publishing  House  for  Metals  and  Science, 
Moscow,  1958. 

93.  S.A.  Saltykov,  Stereology,  H.  Elias  (ed.).  Proceedings  of  the  2nd  Int.  Conf.  Stereology,  Springer- 
Verlag,  NY,  pp  163,  1967. 

94.  T.  N.  Rouns,  J.M.  Fridy,  K.B.  Lippcrt,  0.  Richmond,  Simulation  and  Theory  of  Evolving 
Microstructures ,  M.P.  Anderson  and  A.D.  Rollet  (eds.).  The  Minerals  and  Materials  Society,  Warrendale, 
PA,  pp.  269,  1990. 

95.  Asaro,  R.J.  and  Rice,  J.R.  Strain  localization  in  ductile  single  crystals.  J.  Mech.  Phys.  Solids  25 
(1977),  309-338. 

96.  Anand,  L.  and  Kalidindi,  S.R.,  1995.  The  process  of  shear  band  formation  in  plane  strain 
compression  of  FCC  metals:  effects  of  crystallographic  texture.  Mech.  Mater,  17, 223-243. 

97.  Balasubramanian,  S.  and  Anand,  L,,  2002.  Elasto-viscoplastic  constitutive  equations  for 
polycrystalline  fee  materials  at  low  homologous  temperatures.  J.  Mech.  Phys.  Solids  50,  101  -  126. 

98.  Frost,  H.J.  and  Ashby,  M.F.,  1982.  Deformation- Mechanism- Map,  the  plasticity  and  creep  of 
metals  and  ceramics.  Pergamon  Press,  Oxford. 

99.  S.  Nemat-Nasser,  N.  Luqun  and  T.  Okinaka,  A  constitutive  model  for  fee  crystals  with 
application  to  polycrystalline  OF HC  copper,  Mech.  Materials,  Vol.  30,  pp.  325-341,  1998. 


30 


100.  D.  Pierce,  R.J.  Asaro  and  A.  Needleman,  An  analysis  of  nonuniform  and  localized  deformation  in 
ductile  single  crystals,  Acta  Metallurgica,  Vol.  30,  pp.  1087-1 1 19,  1982. 

101.  Armstrong,  P.J.  and  Frederick,  C.O.,  1966.  A  Mathematical  representation  of  the  multiaxial 
Bauschinger  effect.  CEGB  Report  RD/B/N731,  Berkley  Nuclear  Laboratories. 

102.  E.  Kroner,  Dislocations  and  continuum  mechanics,  Applied  Mechanics  Review,  Vol.  15,  pp.  599- 
606,  1962. 

103.  Needleman,  A  continuum  model  for  void  nucleation  by  inclusion  debonding,  J.  Appl.  Mech.,  Vol. 
54,  pp  525-531,  1987. 

104.  Needleman,  An  Analysis  of  Tensile  Decohesion  along  an  Interface,  J.  Mech.  Phys.  Solids ,  Vol. 
38,  pp  289-324,  1990. 

105.  V.  Tvergaard,  Effect  of  fiber  debonding  in  a  whisker-reinforced  metal,  Mat.  Sci.  Engin.,  Vol. 
A125,pp  203-213,  1990. 

106.  K.  L.  Roe  and  T.  Siegmund,  An  irreversible  cohesive  zone  model  for  interface  fatigue  crack 
simulation,  Engineering  Fracture  Mechanics,  Vol.  70,  pp.  209-232,  2003. 


31 


